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MICROSCOPE GENERANT UNE REPRESENTATION TRIDIMENSIONNELLE D'UN OB JET ET IMAGES 
GENEREES PAR CE MICROSCOPE. 

(57) L'invention concerne un microscope qui calcule une 
representation tridimensionnelle complexe de I'objet a partir 
d'une serie d'enregistrements de I'onde lumineuse issue de 
I'objet et de I'onde eclairant I'objet, I'onde eclairant I'objet 
variant entre chaque enregistrement. 

Le microscope comprend une partie optique permettant 
la generation de figures ^interferences entre une onde lu- 
mineuse de reference et une onde lumineuse diffractee par 
I'objet observe, des capteurs permettant de numeriser ces 
figures d'interferences, des actionneurs permettant d'agir 
sur le systeme optique, et un ordinateur recevant Ies figures 
d'interference numerisees, controiant Ies actionneurs, et 
dote de memoire et de moyens de calcut lui permettant de 
calculer des images tridimensionnelles a partir des figures 
d'interferences. La partie optique permet I'enregistrement 
sur une surface de capteurs (118) de I'onde issue de I'objet 
(112) par mesure des figures d'interferences produites par 
le faisceau de reference Fret I'onde issue de I'objet. Le mi- 
roir (109) permet de faire varier la direction du faisceau 
d'eclairage parallele de I'objet. 

.Application: microscopie en biologie et en metrotogie. 
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Microscopc ecnerant une representa tion tridimcnsionncllc 
dun obiet et images ncncrces par ce microscope. 

Domaine technique: 

5 La prcscntc invention conccrne un microscope optiquc generant unc representation en trois 

dimensions de 1'objct observe. 

Technique antcricure: 

Les microscopes couranls ferment par un precede optiquc une image bidimcnsionncllc 
10 corrcspondant a une coupe aggrandic de 1'objct observe. Cette image peul 1c cas echcanl etre cnrcgistrcc 
par une camera video aftn de pouvoir etre rcstituee ultcrieuremcnt. 

II est possible de genercr une image tridimcnsionnellc en utilisant un de ccs microscopes et en 
faisant varicr le reglage de focalisation. A chaque rcglage correspond un plan de coupe different, et unc 
image en trois dimensions peut ctre reconstiluee a paitir de ccs plans de coupe. Certains microscopes 
1 5 munis d'un dispositif de focalisation motorise et d'un logicicl appropric effectuent cette operation 
autbmatiquemenl. De tels microscopes sont decrits par exemple dans l'article « 4-D imaging software 
observe living cells », par Charles Thomas et John White, paru dans I'exemplaire de decembrc 1996 du 
journal « Scientific Computing World » publie par IOP publishing Ltd a Bristol, Grande-Bretagne. 
II existe egalement des microscopes confocaux, dans lesquels l'eclairage est focalise sur un 
20 point et l'image tridimensionnelle est gendree en balayant tous les points de l'objet. Dc tels microscopes 
sont decrits par exemple dans Particle « Under a Microscope: Confocal Microscopy Casts New Light on 
the Dynamics of Life », par W.Conard Holton, paru en fevrier 1995 dans la revue « Photonics Spectra » 
publico par Optiral Publications Company a Pittsfield. Massachussets, Etats-Unis. 

Par aiUeurs, il existe des microscopes interfcrometriques utilisant deux ondes distinctes dont la 
25 superposition foumit unc image. Sur ces apparcils, l'image est loujours observee directcmcnl dans un 
plan qui est l'image gcometriquc dc l'objet traverse par au moins unc des deux ondes. Ccs dispositifs 
produisent des images bidimensionnelles scmblables a ccllcs produitcs par un microscope ordinaire, a 
ccci pros que les variations d'eclaircment dc l'image ne sont pas caractdristiqucs de l'absorptivite dc 
1'cchanlillon mais des interferences cntre les deux ondes dans un plan image dc l'objet. Ccs dispositifs 
30 nc penncttent pas d'obtcnir d'images cn trois dimensions, lis sont decrits par exemple au chapitre III dc 
l ouvragc « Progress in Microscopy » par M. Francon. paru cn 1961 chc/. Pergamon Press, Oxford, 
Grandc-Brctagnc 

Lc seul moyen permcttani l enregistrcment rapidc d'unc image microscopiquc 
tridimcnsionncllc est la imcroscopic holographique. Mais cette melhodc ne permct pas cn general 
35 1' enrcgistrement direct dc reformation sous forme numeriquc: l'hologramme doit etre enrcgistre sur 
un papier special puis l'image doit ctre reconstiluee optiqucment et observee par un microscope 
classique. En effet le pas d'analyse pour un enrcgistrement holographique doit ctre inferieur a la 
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longucur d'onde. cc qui est liors de portde dcs systemcs dc conversion optique-numenquc actucls. Cettc 
mcthodc est d&ritc par excmple dans Particle intitule « Holographic Microscopy ». par M.Plula, paru 
en 1987 dans Ic volume 10 dc la scne « Advances in Optical and Electron Microscopy ». chc/. 
Academic Press. Londrcs. 

5 Dcs auteurs ont realise des microscopes dans lcsquels unc configuration adaptee dc londe dc 

reference el dc l onde diffuscc par l'echantillon pcrmet un enregistrcment numeriquc. ct dans lcsquels 
la reconstruction holographiquc optiquc est simulec numcriquemenl. Dcs systcmes dc cc type sonl 
decrits dans les articles suivants: 

-Karpov V B: « study of Biological Samples with a Laser Fourier Holographic Microscope >>. Laser 
10 Physicsvoi.4 no3. mai 1994-juin 1994, pages618-623 

-Bianco B ct Al: « Computer simulation of 3D Imaging in Holography : A preliminary study for 
automated Holographic Microscopy ». proceedings of ISMIII '82: first IEEE computer society 
International Symposium on Medical Images and Icons. 24-27 juillet 1984. Arlington.US,pages 232- 
235. . 

1 5 Un microscope cxistc dans lequcl la valeur complcxc cxacte dc l onde lumineuse rccue sur unc 

surface donnee est enregistrcc. dans lequcl la valcur de l'ondc lumineuse dans l'cchantillon est calculcc 
a partir de cet enregistremcnt. ct dans Icquel plusicurs eclairages differents peuvent ctre utilises 
successivement pour supprimcr l'effet dc speckle et obtenir par moyennage une representation 
tridimensionnelle dc l'objet observe. Cependant, ce microscope ne peut obtenir une image tenant 

20 compte dc la phase dc 1'onde que lorsque une seule onde d'eclairage plane est utilisee. La phase dc 
ronde lumineuse dans l'objet. calculee par transformec de Fourier inverse, n est alors pas liee 
directement a l'indice local de l'echantillon mais plulot a son integrate le long du trajct de Tondc 
d'dclairage. Dc plus, la definition dans la direction du faisceau est alors mauvaise. Lorsque ce 
microscope utilise plusieurs eclairages distincts, il les combine par moyennage cn intensitc. ce qui 

25 supprimc reformation de phase. L imagc tridimensionnelle obtenuc est alors cn nombrcs reels et nc 
pcrmet de caracteriser que l'absorptivite de l'cchantillon. Cc microscope est decrit dans lc brevet 
PCT/FR97/01695 

Expose dc l'invcntion; 

30 Une code ctant connue par sa valeur complcxe en chaque point d'une surface de reception, il 

est possible dc genercr sa representation complexe dans une zone observee par un calcul simulant le 
rctour inverse dc la lumicre. La representation ainsi obtenuc est unc fonction complexe dcs coordonnccs 
spatiales. Unc representation frequcnticlle de l onde peut ctre obtenuc par transformation dc Fourier dc 
ccttc representation spatiale. En representation frcquenticllc, les coordonnees d'un point rcpresentcnt la 

35 frequence spatiale (qui est done un vecteur). Unc onde plane a une representation frequcntielle reduitc a 
un scul point dont le vecteur frequence associe a pour normc l'inversc dc la longueur d'onde et pour 
direction la direction de propagation dc l'onde. La representation frequcntielle d'une onde 
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monocluomatiquc ne componc que des frequences de nonnc constantc cgale a I inverse dc la longueur 
d onde ct elle est done bidiincnsionncllc. c^est-a-dirc reduitc a une portion de sphere centrcc sur la 
frequence nullc el dont 1c rayon est Tinvcrse de la longueur d'onde: L/onde issue dc r objet ou Tondc 
d'cclairage pcuvent etrc caracterisecs par leurs representations frequent idles ou spatiaies. Scion le 
5 mcmc principc. toutc fonction complcxc dc la position definic dans la zone observee pcut etrc 
considcree commc une representation spatialeet pcut cgalemcnt ctre convcrtie en representation 
frcqucntielle. On utilisera lc tcrme representation' pour designer une telle fonction indcpendainmcnt 
du domaine spatial ou frcqucntiel dans lequel elle pcut etre cxprimce. 

Le microscope decrit dans lc document PCT7FR97/01695. dans lc mode ou il utilise unc scric 
10 d'ondes d'eclairage distinctcs, genere unc representation tridimcnsionnellc dc 1 objet observe qui nc 
prend en compte que Tabsorptivite locale de Techantillon. L'objcctif de la prcscnte invention est dc 
pennettxe la generation d'unc representation tridimcnsionnellc prcnant egalcmcnt en compie des 
variations dc Tindice de refraction. Pour cela, le microscope doit prendre en compte les differences dc 
phase generics pour chaque eclairage par la presence de Techantillon. Afin dc pouvoir prendre en 
1 5 compte ccs differences de phase, ralgorithme de calcul de la representation tridimensionnclle doit 
disposer d'une representation de chaque onde issue dc robjet ct d'une representation de chaque ondc 
d^eclairagc correspondantc, qui conslituc une reference de phase. L'objet de Invention est un 
microscope qui genere par calcul numeriquc une representation tridimcnsionnellc dc robjet observe a 
partir de plusieurs enregistrcments de Tondc diffractee par robjet et re?ue sur unc surface de reception, 
20 ces enregistrcments different cntre eux par i'onde d'eclairage utilises II differe du microscope decrit 
dans lc document PCT/FR97/01695 en cc qu il comportc des moyens pour determiner une 
representation de chaque onde d'eclairage dans robjel observe. 

Dans la mesurc ou on etudie les proprietes locales de Techantillon observe, on s'interessc 
cssentiellcmcnt au rapport dc Tonde diffractee a londc d'eclairage. dont la phase correspond a unc 
25 difference dc phase induitc par la presence dc robjet el dont r intensite est liee a V intensite absorbec par 
ledit objet (en l'abscence d'objet. ce rapport est constant ct cgai a 1). Scion une version avantagcusc dc 
Invention, ralgorithme de calcul de la representation tridimensionnclle de robjet comporte unc ctape 
au cours de laqucllc il calculc une representation normaliscc de chaque onde issue dc robjet par rapport 
a londc d'eclairage. Une representation normalised est telle quen representation spatialc sa valcur en 
30 tout point soil lc rapport dc Tonde issue dc robjet a londe d'eclairage correspondante. Elle n'est pas en 
cllc-mcmc la representation d'une onde ct sa representation frequenticllc iVest pas toujours 
bidimcnsionnellc. 

Selon unc version avantageuse dc r invention, unc valeur complcxc caractcristique est calcuice 
cn chaque point dc robjet en moyennant les vaicurs des representations normalises obtcnues en cc 
35 point pour chaque eclairage. Cc moyennage peut ctre effectuc, par exemplc, sur renscmblc des 

eclairages produisant au point considere une intensite supericurc a une valcur mirumalc choisie pour 
limiter lc bruit. L'enscmble des valeurs ainsi obtenucs consume la representation tridimcnsionnelle dc 
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l'objet. Ccl.c representation tridimcnsionnclJc prcnd cn compte les differences dc phase induites par la 
presence dc Pobjet et done les variations d'indicc dans l'cchantillon. Dc plus, lorsqu'un noinbre 
suffisant d'cclairage dislincts. choisis dc manicrc appropricc. est utilise, la valcur dc cette 
representation iridimensionnclle en un point tend vers unc liniite dont la partie rccllc depend 
5 esscnUcllemcnl dc labsorptivite locale dc l'cchantillon et donl la partie imaginairc depend 
csscntiellemcnl de l'indicc local dc l'cchantillon. 

Scion une version avantagcusc de ['invention. Ics ondes d'eclairage sont des ondes planes 
different cntrc elles par leur direction, ^utilisation d'ondes d'cclairage planes simplifie l'enscmbic du 
systeme. 

1 o Dans le cas des ondes d'cclairage planes, et lorsquc 1 echantillon utilise est faiblemcnt 

diffractant. la representation frcqucntielle dc londe issue dc l'objet comprend l'onde d'eclairage 
attcnucc. dont la representation est ponctucllc (e'est une frequence pure) ct dont 1'intcnsite est clcvcc, et 
i'ondc diffractcc. dont la representation est rcpartie sur un large domainc de frequences. En premiere 
approximation, la representation frequentielle de londc d'cclairage est done limitcc a un point dont les 

15 coordonnees sont celles du point d'intensite maximalc dans la representation de I'ondc issue.de l'objet. 
la valcur de la representation frequentielle dc 1'onde d'eclairage en ce point etant egalc a la valeur dc la 
representation frcqucntielle dc l'onde issue dc l'objet en ce point. Selon unc version avantageuse dc 
1'invention, la representation frequentielle dc l'onde d'eclairage est done obtenue en extrayant le point 
d'intensite maximale de la representation frequentielle de I'ondc issue de Tobjet. 

20 Dans le cas des ondes d'eclairage planes, et selon une version avantageuse de l'invention, la 

normalisation est eflectuee dans le domainc frcqucntiel en divisant la representation frequentielle de 
l'onde issue dc l'objet par la valcur de l'onde d'cclairage sur la frequence unique de cette ondc, puis en 
translatant la representation frcqucntielle de l'onde issue de l'objet de manicrc a ramener a zero la 
frequence unique dc l'onde d'cclairage. Dans cc cas. la representation normalises de l ondc issue de 

25 l'objet restc bidimcnsionnelle, mais n'est plus ccntrcc sur la frequence nulla 

Dans lc cas des ondes d'cclairage planes. l onde d'cclairage est d'intensite constantc dans 
l'objet. Suivant lc principe indiquc plus haul, lc moyennagc est done effectuc cn chaquc point sur la 
totalitc des ondes d'cclairage. Cctlc operation dc moyennage est lineaire. EUc pcut done ctrc effectuee 
indilTcremmcnt dans lc domaine frequcnticl ou dans lc domaine spatial. Scion unc version avantagcusc 

30 de Tinvcntion. cc moyennagc est effectuc dans lc domainc frcquenticl. ou la representation normaiiscc 
obtenue a partir d un eclairagc donne est bidimcnsionnelle ct non tridimcnsionnelle. cc qui limitc lc 
volume de calculs a effectucr Unc representation frequentielle tridimcnsionnelle de Tobjet est alors 
obtenue et la representation spatialc iridimensionnclle peut cn ctre deduitc par transformation de 
Fourier iridimensionnclle inverse. 

35 La presence dc 1 'objet se traduit, du fait des variations d'indicc et d'absorptivite, par la 

generauon locale d une ondc diffusee qui sc superpose a l'onde d'eclairage. Le rapport de cette ondc 
diffusce localement a l'onde d'cclairage est indcpendanl dc l'cclairage el caractcrise lc point considcrc. 
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La partic recite dc cc rapport caracterisc I'absorphvitc de rcchantillon ct sa panic imaginairc 
caracterisc son indicc. II esl particuliercmcnt intcrcssant dc gencrer cc rapport comme representation 
tridimensionnelle dc robjet. Or la representation gencrce par lc moyennagc precedemment indique nc 
correspond pas a cc rapport. Par cxetnplc. V image gencree par un objet absorbant sur un fond 
5 transparent sera en Mtoir sur blanc* (ondc issue de 1'objct absorbant faiblc. onde intense issue du fond 
transparent fond), alors que lc rapport recherche sera en blanc sur noir (ondc diffuscc nulle sur les 
parties transparent ct elevee sur les parties absorbantcs) . Pour obtcnir ic rapport recherche, ii est 
necessairc de remplaccr le moyennage sur Icnsemblc des representations frcqucntielics normalisces des 
ondes issues de robjet par un moyennagc limiie en chaque point a celles dc ccs representations qui 
10 atteigncnt le point considere. c cst-a-dirc qui ont unc valcur non nulle en cc point. Unc valcur nulle esl 
affectcc aux points non atteints. cc qui est rclativemcnt arbilrairc mais minimise les errcurs qui en 
dccoulent. Ccttc limitation du moyennagc constituc une version avantageusc de V invention. 

Les variations dc direction dc I'onde dc reference, caracterisee par son vectcur frequence, 
pcuvent sc faire de divcrscs maniercs. Le maximum de precision est obtcnu lorsquc toutcs les directions 
15 possibles sont utilisees. Ccpcndant il sufTit de faire varicr le vectcur frequence dc sorte que son 
extrcmite decrivc un ensemble a une dimension pour obtcnir une representation frcquenticllc 
tridimensionncllc. Selon unc version avantageuse dc r invention, son cxtremitc decrit un ccrclc centre 
sur L'axe optique du microscope et dc rayon suffisant. 

Divcrses methodes permettent dc faire varier la direction dc londe plane d'eclairagc. La 
20 solution la plus immediate est d'utiliser un jeu de miroirs plans mobiles permettant de redirigcr un 
faisceau laser parallele d'cclairage vers rcchantillon sous les angles voulus. Mais eclte solution 
ncccssite lutilisation d'un systeme mecaniquc complexc articulant les iniroirs entrc cux. Une solution 
plus simple est d'utiliser un scul miroir renvoyant le faisceau vers un condenseur sous un angle 
variable. Lc condenseur. formant dans le plan observe dc Tobjet une image dudit miroir, redirigc les 
25 ondes d^cclairage vers robjet, sous des angles eleves pouvant allcr jusqu'a son ouvcrture maximalc. 
Ccpendanl. lc condenseur nc transfonnc pas un faisceau parallele en un faisceau parallele. Afin que le 
faisceau soit parallele en sortie du condenseur. il doit etre focalisc en entree sur un point dc la surface 
focale objet dudit condenseur. Selon unc version avantageusc dc l'invention. lc systeme permettant la 
generation des ondes d'cclairage parallclcs comportc un miroir mobile renvoyant le faisceau d^cclairage 
30 vers un condenseur qui lc rcdirige vers la zone obscrvee dc Fobjcl. ct des elements optkres (lcntilles) 
modifiant le faisceau d'cclairagc avanl son entree dans lc condenseur. de sorte que son print dc 
focalisation soil sur la surface focale objet du condenseur. 

L' ouvcrture du faisceau en entree du condenseur determine alors la largcur du faisceau 
parallele sonant du condenseur pour illuminer rcchantillon. Ccitc largeur doit corresponds a la /.one 
35 observee dc rcchantillon. Trop elevee, clle genere du bruit supplemental re du aux interactions de londe 
d'cclairage avee les points dc robjet sc situant hors dc la zone d' observation. Insuirisammcnt elevee. 
elle ne permct pas d^claircr convenablcmcnt la zone observee. Scion une version avantageusc de 
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l'invention, un diaphragme permetiant de regler I'ouverturc dudit faisceau d'eclairage est inscre dans le 
systeme optique faisant converger l'onde d'eclairage vers un point de la surface focaie objel du 
condenseur. 

Selon une version avantageuse de Invention, le microscope est construit de sorte qu'une onde 
5 plane ait une image ponctuelle sur la surface de reception. La valeur de l'onde en un point de la surface 
de reception est alors egale a la valeur de la representation frequentielle de l'onde en un point 
correspondant. 

Selon une version avantageuse de l'invention, une image intermediate de I'echantillon est 
formee et filtree spatialcment par un diaphragme, ce qui supprime les problemes de replicment de 

10 spectre dus a la partie du trajet des faisceaux d'eclairage qui est a l'extericur de la zone observee. 

L/ acquisition d'un enregistrement de l'onde diffractee par I'objet se fait par enregistremcnt 
successif de plusieurs figures d'interferences entre 1 'onde diffractee par 1 'objet et unc onde de reference, 
la phase de l'onde de reference variant entre chaque acquisition et l'onde diffractee etanl calculee point 
par point a partir de ces figures d' interferences. Plus 1'intensile* des figures d' interference est eleven 

15 moins le resultat est brake. Mais lorsque les figures d'interferences deviennent trop intenses en un 
point, le capteur sature localement et le resultat de 1'acquisition est localement invalide. Pour que le 
rdsultat de l'acquisiuon soit localement considere comme valide. le capteur doit etre non sature au point 
considere et en ses voisins immediate, et ce pour I'ensemble des figures d'interferences utilisees. L;onde 
d'eciairage concentree en un point peut etre tres superieurc cn intensite a l'onde diffuses repartie sur 

20 toute la zone de frequences, et lorsque le capteur est non-sature pour l'onde d'eciairage le bruit sur 
l'onde diffusee peut etre considerable. Selon une version avantageuse de l'invention. un dispositif 
d'attenuation contrdlee du faisceau. ayant un ou plusieurs niveaux d' attenuation, est introduit pour 
resoudre ce probleme. Le dispositif d'attenuation permet d'obtcnir successivement plusieurs 
enregistrements different par 1'intensite de l'onde d'eciairage. Une valeur moins bruitee de l'onde 

25 diffractee est alors obtenue en combinant ces enregistrements. La valeur definitive de l'onde diffractee 
est calculee en chaque point a partir de 1' enregistrement pour lequel 1'intensite de l'onde re?uc au point 
considere est la plus elevee. mais pour lequel le capteur reste non-sature au point considere et en ses 
voisins immddiats pour l'ensemble des figures d'interferences permcttant d'obtenir ledit 
enregistrement. 

30 calculs de representation frequentielle tridimcnsionnelle mettent en oeuvrc des quantitcs 

importantes de donnees. Ces donnees etant normalemenl accedees dans un ordre aleatoirc, elles ne 
peuvent pas etre stockees pendant les calculs sur un support a acces sequenlicl comme un disque dur el 
doivent etre stockees sur un support a acces aleatoirc comme une m^moire interne d'ordinateur (RAM). 
Selon unc version avantageuse de I nvention, l algorithme de calcul est modifie de maniere a traitcr les 

35 donnees bloc par bloc, un bloc correspondant a une quantite importante de donnees qui peuvent alors 
etre stockees sequentiellement sur un support a acces scquentiel et chargees en memoirc centrale 
- uniquement pendant le temps de traitemeni dudit bloc. L'algorithme modifie effectue dans un espacc 
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tridimensionnel des iraiiements plan horizontal par plan horizontal, chaque plan horizontal etant stocke 
sur le support a acces sequentiel cn un scul bloc. Afin de pouvoir effectuer egalement des traitements 
suivant la dimension vertical* ralgorithme integrc des phases d'echangc d'axes qui permeltent de 
ramener provisoirement I'axe vertical dans un plan horizontal. La procedure d'echange d'axes opere 
bloc par bloc, les blocs ayant generalemcnt des dimensions egales ou proches selon les deux axes a 
echanger et ayant pour taille en octets la taille maximale qui puisse etre stockee en memoire ccntrale du 
systeme (memoire RAM a acces aleatoire). 

Le microscope decrit genere des representations tridimensionnelles sous la forme de tableaux a 
trois dimensions de nombres complexes stockes sur un support de donnees quel qu'il soil. Selon une 
version avantageuse de l'invention, il g<mere a partir de ces tableaux des images bidimensionnclles de 
Tobjet observe, qui peuvent par exemple etre des coupes ou des projections rcpresentant par exemplc la 
partie reelle ou imaginaire de la representation tridimensionnelle, et qui peuvent etre imprimees ou 
visualisees sur un ecran d'ordinateur. II peut egalement generer des images stereoscopiques constituees 
de deux projections selon des directions difierentes. 



Description sommaire des dessins: 

La Fig. 1 est un schema d'ensemble de 1'optique du microscope. Les Fig.2 et 3 represented le 
detail du positionneur angulaire (1 10) deja represente sur la Fig. 1. La Fig.4 est un schema d'ensemble 
du support mdcanique vertical et antivibratoire du microscope. La Fig.5 represente le detail d'un 
tendeur de la Fig.4. La Fig.6 represente un exemple de dimensionnement de la partie optique. Les Fig. 7 
k 9 sont des representations graphiques servant de support a 1'explication du principc de 
fonctionncment du microscope. La Fig.10 represente l'algorithme d'un programme permettanl de regler 
les tensions de commande du piezoelectrique (122). La Fig.l 1 represente ralgorithme detaille d'une 
procedure de prise d'images utilisee dans le programme precedent. La Fig. 12 represente ralgorithme 
25 d'un programme permettant de regler 1'attenuateur d'intensite constitue du polariseur ( 105) et du 
rotatcur de polarisation ( 104) et d'obtenir ses caracteristiques. La Fig. 1 3 represente I 'algorithme 
detaille d'une procedure de prise d'images utilisee dans le programme precedent et dans les 
programmes d'acquisition d'image 2D ou 3D. La Fig. 14 represente ralgorithme d'un programme de 
focalisation permettant d'obtenir une image 2D et de focaliser l'objectif de microscope (1 13). La Fig. 15 
30 represente l'algorithme d'un programme permettanl de regler la position du condenseur ( 1 1 1 ) et 

d'obtenir certains parametrcs du systeme. La Fig. 16 represente symboliquemenl une image affichec a 
1'ecran par le precedent programme. La fig. 17 represente graphiquemenl un calcul eficctue par cc 
programme. La Fig. 18 represente l'algorithme du programme d'acquisition d'image tridimensionnelle. 
La Fig. 19 represente l'algorithme d'un programme de calcul generant a partir des resultals de 
3 5 l'acquisition une representation tridimensionnelle de 1' objet. La Fig.20 represente schematiquemcnl une 
operation effectuee par la premiere partie de ce programme. La Fig.21 represente le detail de 
l'algorithme de cette premiere partie. La Fig.22 represente schematiquement une operation effectuee 
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par la deuxieme partie de ce programme. La Fig.23 represent l'algorithme d'unc troisiemc panic dc ce 
programme. La Fig.24 representc l'algorithme d'unc dernierc partie de ce programme. 

Description d'un mode de realisat ion uarliciilier. 
5 Scion le mode de realisation prefere. les caracteristiques du microscope sonl telles qu'exposecs 

dans la suite du texte jusqu'au titrc « variantes » non inclus. 

caracteristiques matcrielles: 

Un faisceau laser dc longueur d'onde 633 nm polarise dans le sens orthogonal a la figure est 

10 emis par le laser helium-neon (100) et traverse le filtre (101) constitue de zero, un ou plusieurs filtres 
empiles, en verre teinte Schott. II est ensuite separe en un faisceau d'eclairage Fe et un faisceau dc 
reference Fr par le miroir semi-transparent (102). Lc faisceau d'eclairage traverse ensuite un filtre (103) 
du meme type que ( 101), puis un rotateur de polarisation (104) a base de cristaux liquides 
ferroelectriques et commercialise par la soci&e Displaytech (Euus-Unis), el un polariseur (105). Lc 

15 faisceau traverse ensuite un achromat (106) puis un diaphragme (107) et un achromat (108). 11 est 

reflcchi par un miroir (109) fixe sur un positionneur angulaire (1 10) commande par deux moteurs pas a 
pas. 11 traverse ensuite un condenseur ( 1 1 1) d'ouverture 1.2 et atteint l'objet (1 12). L'objet (1 12) est un 
echantillon place entrc lame et lamella dont l'absorptivite et les variations d'indice sont relativement 
faibles. et dont on veut obtenir une image tridimensionnelle. Les lentilles (106) et (108) doivenl etre 

20 telles que le faisceau laser focalise dans le plan focal objet du condenseur ( 1 1 1 ) et ait une ouverture 
suffisante a son arrivee dans ce plan focal. Le faisceau est done a nouveau parallele lorsqu'il atteint 
l'objet (112). Le plan de l'objet doit etre horizontal pour que l'huile optique necessaire pour utiliser 
l'objectif et le condenseur a immersion nc coule pas. L'axe optique est done vertical et les elements de 
la Fig. 1 sont fixes sur une plaque verticalc (300) representee sur la Fig.4. 

25 Le positionneur angulaire (1 10) est detaille sur les Fig.2 et 3. La Fig 2 rcpresente unc plaque 

sur laquelle le miroir (109) est colle. en vuc de dessous. La Fig.3 represente 1'enscmble du dispositif vu 
de cote. La plaque (200) est posee sur des plots de contact mobiles (205) et (212) dont les points dc 
contact avec la plaque sont en (203) et (204), et sur un plot dc contact fixe dont le point de contact est 
en (201), la plaque (200) etant legcrcment creusce en ce point. Elle est maintcnue par un ressort (210) 

30 fixe a la plaque en (202) et a un plot fixe (2 Si). Le plot de contact mobile (205) se deplacc 

vcrticalement en translation et est integre a un actuateur un axe de type classiquc: il est bloquc en 
rotation par un flexible (207) et est entrainc en translation par la rotation de la tigc (206) sonant du 
moteur pas a pas (208) qui est par cxemple un moteur pas a pas 400 pas/tour. L'entrainement sc fait par 
exemple par T intermediate d'un pas de vis. Le plot de contact (212) est dc meme entraine par la 

35 rotation du moteur (213). Le point de contact fixe (201) doit etre dans l'axe optique de l'objectif (1 1 3) 
et le centre du miroir doit etre sur ce point de contact. Les moteurs (208) ct (213) sont fixes sur une 
plaque (209) elle-meme fixte sur la plaque support (300). 
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L'onde issue de l'objet (1 12) traverse robjectif de microscope (113). Cet objectif est un objectif 
plan (qui donne unc image plane d'un plan), a grande ouverture (par exemple 1,25), a immersion, et 
formant unc image aggrandie de l'objet a une distance finie. 

Dans le plan ou 1 'objectif forme normalement 1'imagede I'objet-a observer, on interpose un 
diaphragme (1 14) permcttant un filtrage spatial de 1'imagc. En arriere de ce plan on positionne un 
achromat ( 1 15) dont le plan focal objet doit etre confondu avec le plan focal image de robjectif (113). 
Un second achromat (1 17) dont le plan focal image est dans le plan d'un captcur CCD (1 18) forme 
dans le plan de ce CCD l'image du plan focal image de I'objectif (1 13). Le CCD (1 18) est intcgre a unc 
camera (119) sortant un signal video analogique. 

Le faisceau de reference traverse d'abord un filtre (120) de mcme type que (101) puis est 
reflechi par un miroir (121) monte sur l'extremite mobile d'un translateur piezoelectrique (122). II 
traverse ensuite une lentille (123) qui focalise le faisceau en un point. Le faisceau divergent issu de ce 
point est reflechi partiellement par le miroir semi-reflechissant (1 16), ce qui le superpose au faisceau 
issu de l'objet et permet d'enregistrer leurs interferences sur le CCD (1 18). Le point de focalisation du 
faisceau issu de la lentille (123) doit avoir son image virtuelle apres reflexion sur le miroir semi- 
transparent (116) au centre de l'image du diaphragme (1 14) par 1 'achromat (115). Le translateur 
piezoelectrique (122) permet de moduler la phase du faisceau de reference. 

Les positionneurs manuels utilises dans le systeme ne sont pas representes sur la figure. 
L'objectif de microscope (113) est monte sur un dispositif de focalisation. Le laser (100) est monte sur 
un positionneur deux axes permettant d'en regler la direction. La partie fixe de la pile piezoelectrique 
(122) est montee sur un positionneur deux axes permettant une rotation par rapport a un axe orthogonal 
au plan de la figure et passant par le centre du miroir (121), et par rapport a un second axe situe dans lc 
plan du miroir (121), orthogonal au premier axe et passant par le centre du miroir (12 1 ). Lc condenseur 
est monte sur un positionneur trois axes en translation. La camera (1 19) est montee sur un positionneur 
6 axes: trois axes de translation et trois axes de rotation passant par le centre du CCD (119). La 
position angulaire des miroirs semi-transparents (102) et (1 16) peut etre ajustee manuellcmcnt. La 
lentille (106) peut etre translatee suivant son axe. L'objet (1 12) est fixe sur un positionneur a deux 
dimensions permettant de le deplacer dans un plan horizontal Le diaphragme (114) peut etre deplacc 
dans un plan horizontal. 

L'ensemblc est fixe air la plaque support (300), du cote de la plaque oppose au point de vue de 
la Fig.4. Cette plaque est fixoe a deux plaques triangiJaires (301) et (302) elles-mcmes fixees a unc 
base carrce (303). La plaque (300) est egalcment fixee directement a la base carree (303). Les plaques 
(300)(301)(302)(303) sont en alliage d'aluminium rigide AU4G, par exemple d'epaisseur 20 mm. La 
fixation des plaques peut se faire par des vis et des trous taraudes, et doit etre faite en un nombre de 
points suffisant pour assurer une rigidite parfaite de l'ensemble. Ceci permet dc maintenir le systeme a 
la verticale en assurant une rigidite suffisante. L'ensemble est pose sur un support antivibratoire 
constitue d'une plaque en granit (304) d'epaisseur 30 mm posee sur une chambre a air de camionclte 
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gonflee a faible pression (305) qui amortit ies vibrations cl qui est clle-mernc posee sur une table rigide 
en bois (306). Un cadre rigide en bois (31 1) est fixe en hauteur par rintermediaire de montants (307) 
(308)(309)(3 10) a la table (306). L'ensembie de la construction en bois est renforcc de manierc a etre 
parfaitement rigide. Le haul de la plaque (300) est relie par des tendeurs (3 12)(3 1 3)(3 1 4)(3 1 5) aux 

5 coins du cadre rigide. Chaque tendcur, detaille sur la Fig.5, est constitue d'un ensemble de bracelets 
elastiques (316) tendus cntre deux anneaux (3 18)(3 17), ccs anneaux etant eux-memes fixes a la plaque 
(300) et au cadre (311) par des cordelettes (320)(3 19), r ensemble etant mis sous tension. La chambre a 
air (305) permet devoir pour rensemble suspendu en AU4G de faibles frequences de resonnance pour 
les mouvements de translation, de Tordre de 2 Hz. Les tendeurs (3 12) a (3 15) permettcnt de limitcr le 

10 balancement de rensemble. La frequence de balancement peut etre evaluee simplement en imprimant 
un leger mouvement de balancier a r ensemble et en mesurant le temps necessaire pour avoir, par 
exemple, dix allers-retour. La frequence de balancement se regie en modifiant le nombre de bracelets 
elastiques utilises pour chaque tendeur. Plus il est eleve, plus la frequence dc balancement augmentc. 11 
doit etre regie pour que la frequence de balancement soil du meme ordre que la frequence de resonnance 

1 5 pour les mouvements dc translation, soit environ 2Hz. 

Le montage des divers dements sur la plaque (300), et en particulier des miroirs et miroirs 
semi-transparents, doit etre effectue de maniere a assurer une rigidite maximale de rensemble. Toutes 
les precautions usuelles doivent etre prises de maniere a limiter les vibrations. 

Les precisions suivantes concernent un exemple particulier de dimensionnement pratique du 

20 dispositif Les distances citees sont representees sur la Fig.6. Les dimensions donnees sont 

approximatives et certaines d'entrc clles doivent etre corrigees dans une phase de reglage. La lentille 
(106) est un achromat de distance focale 10 mm. La lentille (108) est un achromat de diametre 24 mm 
et dc distance focale 120 mm. La distance D7 entre (106) et (108) est dc 250 mm. La distance D6 cntre 
la lentille (108) et ic centre du miroir (109) est de 100 mm. La distance D5 entre lc centre du miroir 

25 ( 1 09) ct le plan focal objet du condenseur est de 1 40 mm. Le condenscur est un condenseur fond clai r a 
immersion d'ouverturc L2. L'objectif dc microscope est un objectif plan xlOO d ouvcrture 1,25, a 
distance finie, formant r image a 160 mm du col de l'objectif. de distance focale environ 1,8 mm. La 
distance D4 entre lc col de r objectif et le diaphragme (114) est de 160 mm. La distance D3 cntre le 
diaphragme (1 14) ct r achromat (115) est dc 20 mm. L'achromat (1 15) a unc distance focale de 200 

30 mm ct un diametre dc 30 mm ct sa face la plus bombee est orients vers le miroir semi-transparent 
(116). L'achromat (1 17) a les memcs caracteristiqucs et sa face la plusbombee est cgalement oricntcc 
vers lc miroir (1 16). La distance D2 cntre les deux achromats est dc 85 mm, permcttant d'inscrcr un 
miroir semi-transparent (1 16) de dimensions suffisantes. La distance entre V achromat ( 117) ct lc CCD 
(118) est de 200 mm. La lentille (123) a un diametre dc 4mm et une distance focale de 6mm. La 

35 distance D9 entre cette lentille et l'axe optique est d'environ 70 mm. La distance D8 cntre V achromat 
(115) et le centre du miroir semi-transparent (1 16), situe sur l'axe optique, est d'environ 45 mm. Le 
laser (100) est un laser helium-neon de longueur d'onde dans le vide A v =633 nm polarise dans lc sens 
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orthogonal a la figure, de puissance environ 0,5 mW, de diam&re de faisceau 0,5 mm. Le caplcur CCD 
est un capteur a pixels carres . la surface du pixel etant d'environ 8,5 x 8,5 micrometres, ct la surface 
utile en pixels £tanl de dimension au moins egale a 512x5 12 pixels. La camera sort un signal video 
CCIR ct une horloge pixel, et son temps d'exposition est egal a la moitie de la duree d'un champ, soil 

5 l/50e dc seconde pour une camera CCIR non entrelacee dont les champs durent l/25e de seconde. Ceci 
permel d'avoir un delai entre la fin d'un champ et lc debut de la periode d'exposition suivantc, delai qui 
pent etre utilise pour modifier les conditions d'eclairage sans que la transition influc sur 1'image. Le 
positionneur pidzoeiectriquc (122) est une pile' piezoelectrique precontrainte en forme de cylindrc dont 
le corps est fixe et l'extremit£ sc deplace de 15 micrometres pour une tension appliqucc de 100 Volts. 

10 Le systcme de calcul est par exemple un ordinateur type TC\ dote de cartes d'acquisition ct de 

commande appropriees et dventuellement de moyens de calcul supplements res, fonctionnant par 
exemple sous le systeme d'exploitation Windows 95. La carte d'acquisition du signal video, 
fonctionnant en temps reel, echantillonnc le signal sur 8 bits et acquiert des images de taille hpix x vpix 
ou hpix et vpix sont superieurs a 5 12 et multiples de 4. Les pixels sont echantillonnes suivant Thorloge 

1 5 pixel, done correspondent exactement aux pixels du CCD. Le positionneur piezoelectrique est pilotc 
directemenl par une carte de conversion digitale/analogique sortant un signal compris par exemple 
entre zero et U max . avec par exemple t/ max =10 volts. Une resistance est intcrposee entre la sortie dc la 
carte de conversion et les bornes de I'actionneur piezoelectrique de maniere a limitcr le courant. Sa 
valeur est reglee de sorte que le temps de montee de la tension aux bornes de I'actionneur (122) soit 

20 d'environ 1 ms. Le rotateur de polarisation (104) est equivalent a une lame demi-onde dont l'axe peut 
tourner et a deux positions d'equilibre separecs par un angle de 45 degrcs. S'il est positionne de 
maniere a ce que le faisceau polarise soit parallele a cet axe dans une des positions d'equilibre. il fcra 
tourner le faisceau de 90 degres dans T autre position. Lc rotateur de polarisation est pilote par 
' application d'une tension bipolaire, -5V correspondant a une position d'equilibre el +5V a Tautrc. 

25 Chaque borne du rotateur dc polarisation est reliee a unc sortie 0/5V d'une carte dc sortie numeriquc, et 
ses deux positions sont commandees en appliquant dans un cas 0V a une sortie et 5 V a Tautre, ct cn 
inversant pour l'autre position. Les moteurs pas a pas (208) el (212) sont egalement pilotes depuis 
Tordinateur, via une carte de commande et une dlcctroniquc appropriees. L'ordinateur est dote d'unc 
memoire interne suffisantc (au moins 32 Mo) et d'un disque dur dc taille suffisantc (au moins 4 Go). 

30 

Conventions divcrscs 

Les conventions suivantes seronl utilisees dans la suite de cctte description: 
-La lettrc j representc tantot un indice. tantot lc nombre complcxc imaginairc pur de module 1. 
Dans le cas ou il peut y avoir ambiguitc lc nombre complexey sera note j . 
35 - le signc = symbolise suivant les cas 1'operation d'affectation ou Tegalite 

-F expression a+=b signifie a=a*-b 
-a multiplic par h est ecrit ah, a.b ou a*b 
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-le module d'un nombrc complexe z est note |z|ct son conjugue est note z . 

- La transformec de Fourier discrete, dans sa forme la plus usuelle, transformc un Dirac situe a 
Torigine en une constante et transforme une constante en un point situe a Torigine. La transformec dc 
Fourier discrete uulisec dans r ensemble de ce brevet transforme une constante cn un Dirac situe au 
5 milieu du tableau transforme et transforme un Dirac situe en ce point en une constante. Ceci signifie 
que le 'zero' en frequence ou en position est place au milieu du tableau et non * Torigine du tableau. 
Cette transformee modifice s'obtient a partir de la forme usuelle en effectuant avant et apres la 
transformation une permutation ^indices Un tableau E de dimension /dim est transforme comme suit: 

1 -premiere permutation: E[i]=E\(i+fdim/2)Vofdim) 
10 2- transformation de Fourier usuelle du tableau E 

3- permutation inverse: E[i]=E\(i+fdiml2)%fdim\ 
ou le signe % designe le modulo. 

-La transformee de Fourier bidimensionnelle d'un tableau de lignes et colonnes est obtenue en 
effectuant la transformation de Fourier monodimensionnelle definie ci-dessus sur chaque iigne du 
1 5 tableau, ce qui genere un tableau intermediaire, puis en effectuant cette transformation sur chaque 
colonne du tableau intermediaire pour obtenir le tableau transforme. 

-De meme. ia transformee de Fourier tridimensionnelle consiste a effectuer successivement 
selon chaque axe des transformers monodimensionnelles etendues a r ensemble du tableau. 

20 Principes de fonctionnement 

La valeur de Tonde lumineuse issue de Tobjet sous un eclairage donne en un point du captcur 
CCD est obtenue a partir de Penregistrement de trois figures d' interferences recmes sur le capteur, la 
phase dc Fonde de reference etant decalee de 120 degres entre chacunc de ces figures. 

Si .v est la vibration lumineuse issue de Tobjet et r est la vibration lumineuse constituant Pondc 
25 de reference lors du premier enregistremcnt. les vibrations lumineuscs totales parvenant au capteur lors 

2jz .In 

/"-- -J 

des trois enregistrement successifs sont: s 0 = s + r , jr, =s + re 3 , s 2 =s + re 3 . Les intcnsites 
enrcgistrees successivement sont done: 

kol 2 = N 2 ^-H 2 



|2 I |2 i |2 



+ 



3 j_ c<rf> 3 



sre 5 +sre 

V J 



f 2n 2n\ 

•3« h| 2 =N 2 +H 2 + 



sre 3 +sre 3 



On peut inverscr ces formules et on obtient: 
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La valeur ci-dcssus est la vibration lumineuse provenanl de l'objet scul, la reference de phase 
etant conventionnellement egale a zero pour une vibration en phase avec l'onde de reference. Ce calcul 
pennet done de reconstitucr la valeur complexe de Fonde a partir des enregistrements d'intensite. 
5 Par construction, chaque point du capteur CCD correspond a une frequence pure f c dc Kondc 

provenant de Techantillon. On appelle centre optique du capteur lc point du capteur illumine par un 
rayon entrant dans 1'objectif avec une direction strictement parallcle a Taxe dc symetrie dc ce dernier ct 
on noteC r C\, ses coordonnees en pixels. L'objectif veriftant la condition des sinus, la deviation en 
sortie d'objectif d'un rayon originaire d'un point central de l'objet est proportionnelle au sinus de son 

10 angle d'entrde dans l objectif , qui vaut sin? = . V y , ou [f x ,f y ,f 2 ) sont les coordonnees 



du vecteur frequence spatiale du faisceau en entrde de Tobjectif, de normc -j ou X est la longueur 

d'onde dans le milieu observe, et de direction la direction du faisceau. Le reste du systeme optique etant 
paraxial, la deviation du point d'arrivee du rayon sur le capteur par rapport au centre optique du capteur 
est egalement proportionnelle a cette grandeur, et done les coordonnees {/ - C x J - C y } de ce point par 

15 rapport au centre optique sont proportionnelles a {f x >f y ) Dans le ^ ou le capteur a des pixels carres, 
le vecteur frequence en entree de l'objectif du rayon qui illumine le point du capteur de coordonnees 

(/,/) vautdonc: f c {ij) = -C xtJ - C y ^K 2 -(/ -C x ) 2 -(f-C y ) 2 } 

ou K est une constante a ddterminer. On appellera 4 frequence caracteristique 1 du point ce vecteur 
frequence. 

20 Connaissant la frequence caracteristique de chaque point du capteur et la vibration lumineuse 

issue dc l'objet et recue en ce point, on obtient done, pour un eclairage donnc de Techantillon, la 
representation frequentielle de 1'onde issue de cct echantillon. 

Lorsqu'un echantillon suffisamment fin, peu absorbant et de faiblcs variations d'indice est 
traverse par un faisceau laser parallcle, chaque point de Techantillon est soumis a une vibration 

25 lumineuse Ae jl7 ^ - r ou f e est lc vecteur frequence spatiale du faisceau d'eclairage et r le rayon- 
vecteur au point considere de l'objet. Torigine etant prise au point d'origine virtuel de Ponde de 
reference. L'absorbtivitc et les variations d'indice de r echantillon se traduisent par r apparition d'une 
onde secondaire qui se superpose a l'onde d'eclairage. Un petit volume dVcrec une onde secondaire 

Ae j2/ ^ e r u(r)dV ou u(r) est un coefficient complexe dont la partie reelle est liee a Tabsorptivite 
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localc de l'echantillon et la partie imaginaire a son indice. L'imagc iridimensionnelle de 1'cchantillon 
que ce microscope genere est I'ensemble des vaieurs u{r) en chaque point de rechantillon. 

Lorsque un point de rechantillon, de rayon-vecteur r . emet localement une vibration 
lumineuse s(r) , la vibration lumincuse recue sur le capteur en un point de frequence caracterislique / B 

5 vaut alors s(r)e~ J2 ^ cr . Un petit volume dVte Tobjet. de rayon-vecteur r. eclaire par une onde plane 
de frequence f e . cree done au point du capteur dc frequence caracterislique f c une vibration 
Ae jl7^f e -f c )j- u ( r )dV qui se superpose a la vibration principal. Integree sur I'ensemble de l'objet. 
la vibration recue en un point du capteur vaut done s(P) = j]J Ae j2 ^ fe ~ fc ^ r u(r)dV . Cettc 
vibration est done un element de la transformee de Fourier de la fonction «(r) , correspondant a la 

10 frequence de Fourier f,=f e - fc Le principe de ce microscope est d'enregistrer cette vibration pour un 
ensemble de frequences f e et f c , puis de reconstituer la representation frequentielle de u(r) ct 
finalement w(r) en inversant la transformee de Fourier. 

Les echantillons etudies sont peu absorbants et ont de faibles variations de I'indice de 
refraction. Par consequence l'onde d'eclairage reste tres intense et illumine un point du CCD, que Ton 

15 appellera point d'impact de l'onde d'eclairage. et qui est le point du CCD ou la vibration lumineuse 
recue est la plus elevee. Si ce point a pour coordonnees (Jmaxjmax) alors la frequence de l'onde 
d'eclairage est: 

f e = J_ ( imax - C x . jmax - C y . fe*- (imax - C x ) 2 - [jmax ~C y ) 2 \ 

Dans les equations ci-dessus la reference de phase est la phase de Tonde d'eclairage au point 
20 de rayon-vecteur r=0. Or dans une representation frequentielle bidimensionnelle acquise directemeni 
sur le capteur CCD, la reference de phase csl la phase de Ponde de reference cn ce point. Dc plus, dans 
les equations ci-dessus. on suppose constante Tintensite de l'onde d'eclairage. II est done ncccssaire de 
normaiiser la representation bidimensionnelle en phase et en intensite par rapport a l'onde d'eclairage, 
pour se ramcner a des conditions standard qui sont une phase nulle et une intensite unitairc de cette 
25 onde d'eclairage au point r=0, e'est-a-dire que le point correspondant de la representation frequentielle 
doit avoir une valeur de I. Ceci est fait en divisant I'ensemble de la representation frequentielle 
bidimensionnelle par sa valeur au point dc coordonnees (imaxjmax). 

La figure 7 montre I'ensemble (500) des frequences caracteristiques correspondant aux points 

du capteur. U s'agit d une portion de sphere de rayon ^ , limitee par rouverturc dc robjectif, ccntrcc 

30 sur T axe optique (50 1 ) du systeme. Un exemple de vecteur frequence d'eclairage (502) est superpose a 
cet ensemble. Enfin, rensemble(503) des frequences /, =f c ~f e correspondantes en est deduit. 
Lorsque Ton fait varier les frequences d'eclairage selon un cercle (504) comme indique sur la figure 8. 
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Tcnsemble dcs frequences enregistrces est celui genere par la rotation de la portion de sphere (503) 
autour de Taxc (501) dans un plan horizontal, comme Tindique la vue cn coupe dc la Fig.9. On obtient 
done une representation frdquentielle tridirnensionnelle (et non plus bidimensionnellc comme dans le 
cas ou un seul enregistrement est utilisd). A partir de cette representation tridirnensionnelle dans 
5 Pcspace des frequences, on peut generer la fonction u(r) par transformation de Fourier inverse. 

L'operation de normalisation decrite dans la description de l'invention est constitute, dans ce 
mode de realisation, d'une pari de la normalisation cn phase et intensite par rapport a l'onde 
d'eclairage dans le domaine frequentiel, d'autrc part par la translation de vecteur -f e permettant 
d'obtenir les frequences f,=f c -fe- En effet > ^operation de normalisation en phase et intensite 
10 ramene a 1 la valeur de Tonde d'eclairage au point de rayon-vecteur nul, ct la translation de vecteur 
~f e ddmodule cette onde d'eclairage et en ramene a 1 la valeur en tout point. Ceci equivaut 
effecuvement k diviser en tout point de la representation spatiale l'onde lumineusc reconstruite dans 
I 'echantillon par l'onde d'eclairage au meme point. 

L'operation de moyennage decrite dans la description de l'invention consiste, dans ce mode de 
15 realisation, a calculer la representation frdquentielle tridirnensionnelle comme moyenne des 

representations frequentielles bidimensionnelles normalises. Cependant on n'obtient pas de cette 
maniere la representation frequentielle de u(r). Pour obtenir la representation frequentielle de u(r) il 
faut limiter le moyennage en chaque point £ Tensemble des representations frequentielles atteignant ce 
point, de maniere a ne pas sous-evaluer les valeurs associees aux points qui sont atteints rarement. 
20 Cette methodc pose un probI6me pratique qui est que le point eclaire dircctement par le 

faisccau traversant l'cchantillon est illumine de mantere beaucoup plus intense que les points 
correspondant a Pondc diffractee. La representation tridirnensionnelle generee lors d'une prise d'image 
contient essentiellement dcs frequences proches de celle du faisceau d'eclairage, les autrcs frequences 
etant noyecs dans le bruit. Pour remedier a cet inconvenient, on utilise un dispositif d 1 attenuation 
25 controlee du faisceau. La panic de la representation frequentielle correspondant aux frequences sur 
lesquelles V intensite est elevec est obtenue avee une forte att6nuation et cellc correspondant aux autres 
frequences est obtenue avec une faible attenuation. Les valeurs obtenues sous forte attenuation sont 
ensuitc multipliees par un coefficient complexe caracteristique du decalage de phase et du rapport 
d'ampiitude de Tonde d'eclairage entre les deux positions du dispositifd' attenuation controlee. Cc 
30 dispositif d'atfenuation controlee est consume du rotateur de polarisation (104) et du polariseur (105). 

La representation tridimensionneHc obtenue a partir d images dc taille hpix x vpix correspond 
a des tallies de fichiers importantes. Afin dc limiter la taille des fichiers ct le temps dc calcul, la taille 
des images sera divisee par deux par moyennage lors dc la procedure d'acquisition dc la representation 
tridirnensionnelle. Ceci equivaut a grouper les pixels 4 par 4, un groupc de 4 pixels sur l'image 
35 originate etant equivalent a un pixel effectif utilise pour le calcul. La taille dc Tobjet observable est bien 
sur reduite en consequence. Les valeurs de C x , C y et K sont divisccs par 2 pour tenir compte du 
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nouveau sysleme de coordonnees. Cette limitation de la taiile d'image observes est bien entendu 
facultative. 

Avant de pouvoir utiliser ce microscope, divers reglages dotvent etre eficctues: reglagc des 
positionneurs manuels, obtention de parametres necessaires a la prise d'image ou au calcul, reglagc du 

5 dispositif d'attenuation controlee. L'utilisation du microscope proprcment dite comporte une phase de 
focalisation sur Techantiilon. suivie d'un dernier reglage, puis une phase d^acquisition des 
representations frequenuelles bidimensionnelles normalises, une phase de calcul tridimensionncl et 
une phase de visualisation. Les r6glages des positionneurs manucls permettent de regler conectemcnt la 
trajecloire du faisceau, de s'assurer que T image d'unc onde plane sur la camera soil effectivemenl 

10 ponctuclle que le condenseur forme bien a sa sortie un faisceau parallelc. Le reglagc des tensions de 
commande du piezoelectrique permet devoir des decalages de phase correspondant effectivemenl a 120 
degres. Le reglage de Tattenuateur d'intensite permet d'avoir une attenuation du niveau necessairc, 
correspondant a une amplitude du signal divisee par 8 (intensite divisde par 64), et d'obtenir la valeur 
exacte du coefficient complexe caracteristique de cette attenuation. Le reglage fin de la position du 

15 condenseur permet de s'assurer que Textremite du vecteur frequence de Tonde d'eclairage puisse 
effectivement parcourir le cercle indique plus haut sans que le faisceau sorte de la zone observee. 
L'onde de reference doit etre reglee a un niveau convenable et enregistree. L'onde d'eclairagc doit 
egalement etre reglee a un niveau convenable avant de lancer l'acquisition. 

20 Manipulation des filtres 

Les operations de reglage necessitent des manipulations permanentes des filtres pour ajuster 
Tintensite recue sur le capteur. Ces manipulations ne sont pas systematiquement rappelees dans la suite 
du texte. 

Lc filtre en (120) determine I^intensite dc Tonde de reference. Sa valeur est ddtcrminee dans 
25 une etapc particuliere du reglage. Par la suite, lorsqu'une onde dc reference est necessairc, le filtre ainsi 
determine est inscre. Lorsque Tonde de reference doit etre supprimee, on insere un element opaque en 
(120). 

Le filtre en (103) determine Tintensite de l'onde d'eclairage. Sa valeur depend des operations 
de reglage en cours. Pour la plupart des operations, le filtre est ajuste de manicrc a ce que Tintensite 
30 recuc sur le capteur (118) soit elevde. sans attemdrc la saturation. Pour certaines operations le capteur 
est sature. Pour d'autres Tondc d'cclairage doit etre supprimee, ce qui sc fait par insertion d'un element 
opaque. 

Lc filtre en (101) est utilise uniquement pour les operations de reglage neccssilant un suivi 
visucl du faisceau. repcre par sa tache de diffusion sur un morceau de papier blanc. Le filtre est alors 
35 regte pour que la tache de diffusion soit visible sans etre dangereuse pour Tocil. 

Dans les cas ou seule Tonde d^eclairage est presente, on appellera intensite relative de Timage 
recue sur le capteur CCD le rapport de Tintensite recue sur le capteur a Tintensite en sortie du filtre 
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(103). Dans une operation ou on cherche a maximiscr Pintensite relative recue sur le capteur, on est 
changer de filtre regulierement pour maintenir r intensity a un niveau mesurable par le 



amene a 
capteur. 



5 Programmes (T usage courant 

Certains programmes simples sont utilises frequemment pendant le reglage, sans que cela soit 

rappele: 

- deplacement du miroir (109): ce miroir etant motorise, un programme est necessaire pour cn 
modifier la position. Ce programme demande a Putilisateur un nombre de pas et un numero d^axe 
10 correspondant soit au motcur (213) soit au moteur (208), puis fait effectuer le nombre de pas demande a 
ce moteur. 

-Visualisation de r image recue sur le capteur: un programme permct I'affichage direct sur 
Tecran de Tordinateur de r image recue sur le capteur (1 18). 

-Visualisation de V image et caracteristiques du maximum: ce programme effectue I'affichage 

1 5 direct sur Tecran de l'ordinateur de l'image recue sur le capteur (118). II affiche en plus la valeur 

maximale detectee par le capteur, les coordonnees du point correspondant, et le rapport entre Tintensite 
de ce point et la somme des intensites de ses 8 voisins. Ce programme est utilise pour verifier 
Tapparence d'une image, pour verifier la non-saturation du capteur (valeur maximale inferieurc a 255), 
pour connaitre les coordonnees et la valeur du maximum, pour apprecier le caracterc ponctuel de ce 

20 maximum par observation directe de l'image et par utilisation des valeurs affichees: Tintensite 

(relative) du maximum doit etre aussi elevee que possible, ainsi que le rapport de son intensite a celle dc 
ses voisins. 

reglage de la position du laser (100) et du miroir (121) 
25 Dans un premier temps l'ondc d'eclairage est supprimee et la position du laser (100) est rcglec 

de facon a viser effectivement le centre du miroir (121), ce qui se verific en suivant le trajet du faisccau 
a i aide d^un morceau de papier permettant de le visualiser. La position du miroir (121) est ensuite 
ajustec de maniere a ce que le faisceau de reference traverse eficctivement la lcntille (123) cl parvienne 
a la camera. Le faisceau de reference doit etre centre sur le capteur (118). 

30 

reglage en translation dc la position de la camera et reglage du m iroir (102) 
La position de la camera est rcglec en translation en envoyant directement un faisceau parallele 
sur Tobjectif dc microscope. A cet effet les elements (106)(108)(1 1 1)(105)(104) sont provisoiremenl 
cnleves, Tonde dc reference est supprimee. Tobjectif de microscope est positionne en position a peu pres 
35 focalisee sur Tobjet. L'objet ( 1 12) utilise est une lame transparentc et de l'huile optique est intcrposee 
entre (1 12) et (1 13). La position angulaire de (102) est alors ajustec pour que le faisccau parvienne 
directement au centre du miroir ( 109). La position du miroir (109) est reglee de sortc que le faisccau 
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parallelc cntre dircctement dans l'objectif (1 13) et est aflinee de maniere a maximiser I'intcnsite 
relative du signal rcca sur le capteur CCD. La position du capteur CCD est alors rcglce en translation 
dans la direction de l'axe optique dc maniere a ce que l'image produite soil parfaitement ponctuelle. 
puis est reglee en translation dans les directions orthogonales a l'axe optique de sorte que cc point soil 
5 au centre de la zone utile du capteur. Ce point central est le centre optique ct ses coordonnees cn pixels 
sont les constantes C x , C y . 

reelaee de la position du conde nseur till') 

Us elements (106)(108)(1 1 1) sont remis en place. De l'huile pour microscope est interposec 
10 entre (111) et (112) et entre (1 12) et (113). Un morceau de canon blanc est pose sur le miroir (109) 
pour difiiiser la lumiere. Le diaphragme (107) est ouvert au maximum cl l'ouverture du diaphragme 
( (114) est d'environ 8 mm. L'objectif de microscope est mis en position a peu pres focalisee. La position 

du condenseur (1 1 1) est alors reglee de maniere a obtenir sur le CCD un disque clair de rayon elevc, le 
disque devant restcr relativement homogenc (les bords ct le centre sont toutefois un peu plus eclaires). 

15 

reelaee en rotation de la position de la camera et r e elaee de la position de la lentille (106) 
Le morceau de carton blanc pose sur le miroir ( 109) est supprime. La posiuon du miroir est 
modifiee de maniere a amencr le point eclaire sur le bord du disque precedemment obtenu, en un point 
de coordonnees (x.O) avec xX), l'axe x etant pris seion un axe situe dans le plan de la figure 1 et l'axe y 

20 selon un axe orthogonal a ce plan. La position de la lentille ( 1 06) est alors ajustee selon son axe de 
maniere a avoir sur le capteur CCD 1'image la plus ponctuelle possible. La position de (109) est^alors 
modifiee pour amener le point eclaire sur le bord oppose du disque. La posiuon de ( 1 1 9) est alors 
ajustee en rotation par rapport k l'axe orthogonal au plan de la figure, dc maniere a amdliorer la 
. ponctualite. cl la position dc ( 106) est egalcment ajustee. Une partie de l'ajustement doit se fairc par la 

25 position de (106) ct une partie par la position de (1 19). Le point eclaire est alors ramenc sur le bord 
initial du disque et la mctne procedure est repetee alternativement sur chaquc bord jusqu'a ce que la 
ponctualite soil optimale des deux cotes, ce qui est obtenu tres rapidement en pratique. Le point eclaire 
est alors amene sur le bord haut du disque. en un point dc coordonnees (0,y) avcc.v>0, et la position 
angulaire de la camera est ajustee pour avoir la meillcurc ponctualite possible. Le point eclaire est alors 

30 ramenc vers le bord oppose du disque pour verifier la ponctuuSite et evcntuellement operer unc derniere 
correction. A ce point, la position de la camera est entieremcnt ajustee. 

Determination du nombre de pas oar pixel . 
Le morceau de papier occultant le miroir est alors enleve. Le moteur est deplace suivant un axe dun 
35 nombre de pas connu. La position en pixels du point d'intensite maximale est notee avant ct apres lc 
deplacement. et le nombre de pixels parcourus en est deduit. On calcule alors lc rapport 
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pas_ par_ pixel - 



nombrc de pas 



. On effectue dc meme sur l'autre axe ct on retient le plus petil 



nombre de pixels 



rapport obtenu. 

rfolage du niveau de l'onde de reference 

Si on considere une onde de reference et une onde d'eclairage ayant meme intensitc maximalc 
£k leur arrivee sur le capteur, s'additionnant en amplitude lorsqu'elles sont en phase, la condition de 
non-saturation du capteur est que l'amplitude commune des deux ondes soit la moitie de l'amplitude 
saturant le capteur, ou de maniere equivalente que l'intensite commune des deux ondes soit le quart de 
l'intensite saturant le capteur; Pour regler le niveau de l'onde de reference a cette valeur, l'onde 
d 'eclairage est supprimee et la valeur du filtre ( 120) est ajustee pour obtenir une image dont le niveau 
maximal est d'environ le quart du niveau maximal autorisd par la carte d'acquisition, soit dans le cas 
d'un echantillonnage sur 8 bits du signal video, un niveau d'environ 64. Avant la prise damages, le 
niveau maximal de T onde d'eclairage devra etre regie de la meme fa<?on. 

reelaee des tensions de commande de Tactuateur piezoelectrique 

L'etape suivante consiste a r^gier les tensions de commande de l'actuateur piezoelectrique. 
Ceci est fait en utilisant un programme qui fait varier les tensions de commande et qui est decrit par 
Falgorithme de la figure 10. Avant de lancer le programme, la position du miroir (109) doit etre reglee 
de maniere a ce que le point produit sur le capteur CCD en utilisant pour objet une lame completement 
transparente soit au centre du capteur. Pour utiiiser ce programme, V objet utilise doit etre fortement 
diffusant. On pourra par exemple utilise un morceau de papier blanc trempe dans de la gdlatine puis 
place entrc lame et lamelle. Le papier doit etre suffisamment epais pour arrcter le faisceau direct et 
suffisamment fin pour laisser passer un faisceau diffuse. Le diaphragme (107) est ouvert au maximum 
et le diaphragme (114) est regie pour une ouverture d'environ 0,8 mm. L'onde de reference interfere 
alors sur le capteur CCD avec le disque clair produit par l'onde issue de l'objet. L'intensite de l'onde 
d'eclairage doit etre reglee pour que le capteur soit en limite dc saturation. 

Si la tension maximale appliquee a l'actuateur est Umax, les tensions correspondent aux 

decalages de phase dc --^,0, seront respectivement Umax/2-diff bas, Umax/2, Umax/2 -\ 

diffhaut, ou diffhaut et diffjbas sont choisk pour produire les decalages de phase indiques. Afin 
d'eviter tout effet d'hysteresis, a chaque acquisition la tension est initialisee a 0, les differentes images 
sont acquises dans I'ordre croissant des tensions appliquees a l'actuateur et unc tension finale Umax est 
finalement appliquee, de sorte que le meme cycle est toujours utilise. Ce cycle devra egalement etre 
utilise en phase normale de fonctionnement. Les images prises avec les decalages de phase indiques 
plus haut permettent de calculer une representation frequentielle en appiiquant a chaque pixel P la 
formule 
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S(P) = 



dans laquelle 1'onde de reference a ete remplacee par une constante et t'expression l{P,a) designe 

l'intensite enregistree au point P pour un decalage de phase a. 

Le reglage des tensions de commande de I'actuateur consiste a evaluer diffjaut et diff bas. 

5 Le principe est dobtenir deux representations frequentielles decalees entre elles de j . La deuxiemc 

ft 

image pent etre recalee en phase en la mulUpliant par i* * et I'ecart moyen entre la premiere image et 
cette image recalee peut etre calcule. Cet ccart est minimal lorsque les tensions diff. haul et diff bas 
sont correctement reglees. Les decalages dc phase permettant l'acquisition de la 2eme image 

elemental en frequence sont done + ^ , f . }j + f et correspondent respectivemenl, en 

10 premiere approximation, aux tensions Umax/2-diff_bas+diff_bas/2, Umax/2+(d.ifr_bas+di/f_haut)/2, 
Umax/2^ diff haul +diff haul/2. Le programme de reglage calcule l'ecart moyen pour des series de 
valeurs de diff haul et diff bas et choisit celles qui correspondent a un ccart moyen minimal. Son 
algorithme est detaille sur la figure 10. 

Les etapes essentielles de cet algorithme sont: 

1 5 (600): acquisition des images. La procedure d'acquisition est precisee sur la Fig. 1 1 . Par * les 

images' on entend ici 6 figures d'interferences repues consecutivement par le capteur CCD. Cette 
procedure effectue toujours le meme cycle commencant par la tension 0 et finissant par la tension 
Umax. Le temps d'attente apres Implication d'une tension nulle au pie2»61cctrique permet a celui-ci dc 
se stabiliser. Le temps d'attente avant lc lancement de Facquisition evite d'acquerir une trame qui 

20 aurait cte exposee avant l application des conditions d'exposition voulues. La fin de l cxposition d'unc 
image est signalee sur les cartes d'acquisition par le debut du transfert de l'image correspondante. 
L'utilisation d'un temps d'exposition inferieur au temps de transfert de l'image permet a l'image de ne 
pas etre affectee par les &ats transitoires. La mise du processus en priorite maximale evite la 
perturbation de l'acquisition par d'autres laches systeme sous un systeme d'exploitation multitaches. 

25 Les 6 images sont acquises successivement et en temps r£cl (pas d'image perdue) par la carte 

d'acquisition afin de minimiser le temps pendant lequcl des vibrations peuvent affecter le resultat. 
Chaquc image a pour dimension horizontale hpix et pour dimension verticale vpix. Pendant 
l'acquisition, les images sont transferees automatiquement par la carte d'acquisition dans le tableau qui 
leur est reserve en memoire centrale de 1'ordinateur. A Tissue de la procedure d'acquisition, on dispose 

30 d'un tableau I\a,b,i,j). l'indice a corrcspondant a la difference de phase, l indicc h correspondent a 
l'image (decaiee ou non decalee en phase), les indices / ety etant les coordonnees en pixels et variant 
respectivement de 0 a hpix-\ et de 0 a vpix-l. 
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(601) calcul des deux representations frequentielles. Celles -ci sont stockees dans des tableaux 
de complexes S0li,j\ et SI en appliquant les formules: 
*°M - ^(240,0,/. j] ~ /[MM./] - 42,0,/,;]) + 7^(/[l,0,/,y] - /[2,0,/,y]) 

sA'J] = i(2/[o.i ! /,y]-/[u,/ ! ;]- i[2Xij})+7 ^(/[u,/.yl-/[24,/jD 



(602) Le programme modifie le tableau SI en multiplianl chacun de scs Elements par e 

(603) : le programme calcule la valeur maximale mod_max = 0 ^^ ix _^ S $^\ du inodule 



0£/£i'/?/x-l 



sur 50. 



(604): Le programme calcule F ecart moyen enire les deux tableaux de la maniere suivante: 
Le programme initialise ecort et nombre vaieurs a 0 et parcourt r ensemble des points ij en 
10 testant la condition (605): |S0[/\ j}\ > 0,5 mod max . Chaque fois que cette condition est realise, il 
effectue les operations suivantes (606): 
ecart+ = \SC[i,j]-Sl{ij]\ 
nombre vaieurs + = / 

Lorsque le programme a termine de parcourir les indices (/',/) il divise ecart par nombre _valeurs (607), 
15 ce qui lui donne Fecart moyen. Cet ecart n'integre done que des vaieurs pour lesquelles Fonde 

d'eclairage est sufSsammeni forte, afm d'eviter un resultat trop bruite. II est stocke dans un tableau. 

(608) : le programme passe a la valeur suivante du couple (difj \ haul Jiff _bas) et reitere les 
operations jusqu *k ce que la serie d'ecarts ait ete calculee. 

(609) : le programme applique un filtre passe-bas au tableau obtenu pour limiter le bruit du aux 
20 vibrations du systemc. Le filtre passe-bas utilise est represente dans le domaine frequentiel par une 

marchc d^scalier 1 passant d'une valeur dc 1 en basse frequences a une valeur de 0 en hautes 
frequences, ct sa bande passante est determinec empiriquement pour avoir unc bonne limitation du bruit 
sans trop deformer la courbe. 

(610) : le tableau est represente graphiquement pour cn verifier failure. 

25 (611): la valeur du couple (diff_bas,diff_haut) corrcspondant au minima est affichcc. 

Les series de vaieurs du coup\e(diff_bas t diff_haut) sont determinecs comme suit: 
Dans un premier temps le programme fait varier les vaieurs de diffhaut ct diff_bas en les 

Umax U max 

laissant egales entre elles. Par cxemple elles peuvent varier entre 0 et — ^— par pas de si on 

utilise 12 bits de conversion. II en resulte un tableau de 1024 elements (les ecarts calcules pour chaque 
30 valeur de (diffjiaut.diff_bas) ) que le programme filtre pour eliminer le bruit (609) et qu'il represente 
graphiquement (610). La valeur d'indice du tableau correspondant a Fecart minimal correspond alors a 
la valeur correcte du couple (diff_bas t diffhaut) et ce couple est affichc (61 1). 
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Dans un deuxieme temps, le programme est relance en fixant diffjas a la valeur 
precedemment obtenue et en faisant varier uniquemenl diffjaut. Un nouveau minima est ainsi obtenu, 
qui correspond a unc valeur plus precise de diffjaut. 

Dans un troisieme temps, le programme est relance en fixani diffjtaut a la valeur 
precedemment obtenue et en faisant varier uniquement diffjas, oblenant une valeur plus precise de 
diffjas. 

L operateur pcut reitcrcr ces etapes en faisant varier separernent et alternativement diff has el 
diffjaut, mais le maximum de precision sur ces valeurs est obtenu assez rapidement. 

reglaee du polariseur (105) ct du rotateur de pola risation (104) 

L'etape suivante consiste a regler la position du rotateur de polarisation (104) et du polariseur 
(105). Cc dispositif est destine a rcaliser une attenuation contrdlee du faisceau d'eclairage par 
commande du rotateur de phase, et on appeilera l'ensemble 'commutateur optiqueMl a une position 
fermec correspondant a une intensite faible le traversant et une position ouverte corrcspondant a unc 
intensite plus elcvee. 

L'objet utilise est 1c meme que dans Tetape precedente, le reglage du miroir (109) ct des 
diaphragmes est egalement le meme.L'onde de reference est d'abord supprimee. Dans un premier 
temps, le polariseur est mis en place et regie pour maximiser 1' intensite lc traversant ct recue sur lc 
capteur (1 18). Dans un deuxieme temps, le rotateur de polarisation est mis en place. Une tension 
corrcspondant a un etat arbitraircment defini comme ferme lui est appliquee et il est positionne en 
rotation pour minimiser r intensite traversant r ensemble du commutateur. Dans un troisieme temps, 
Tonde de reference est retablie ct on utilise un programme calculant le rapport des intensites et la 
difference de phase entre les deux positions du rotateur de polarisation, rapport qui sera neccssaire 
pendant les phases d'acquisition d'image pour combiner les ondes correspondant aux etats ouven et 
ferme du commutatcur. L'algorithmc dc ce programme est sur la figure 12. Les dlapes sont les 
suivantes: 

(800) : Le programme acquicrt d'abord les images en utilisant la procedure acquiert images 
decrite Fig. 1 3. Par 'les images 1 on cntend ici 6 figures d' interferences recucs consecutivement par lc 
capteur CCD. Cette procedure effectue toujours lc meme cycle commencant par la tension 0 ct finissanl 
par la tension Umax. Le temps d'attente apres Implication d une tension nulle au piezoelectrique 
perniet a eclui-ci dc se stabiliser. Lc temps d attcntc avant lc lancement dc r acquisition evite d'acqucrir 
une tramc qui aurait 6t6 exposee avant Implication des conditions deposition voulucs. La fin dc 
Imposition d une image est signalec sur les cartes d'acquisition par le debut du transfer! de Timagc 
correspondante. L'utilisation d'un temps deposition inferieur au temps de transfert de Timage permci 
a Timage de ne pas etre affectee par les etats transitoircs. La mise du processus cn priorite maximale 
evite la perturbation de l'acquisition par d'autres taches systeme sous un systeme d' exploitation 
rnultitaches. Les 6 images sont acquises successivement et en temps reel (pas d'imagc perdue) par la 
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carte d' acquisition afin de ininimiser le temps pendant lequel des vibrations peuvent afTectcr le resultat. 
Chaque image a pour dimension horizontal hpix et pour dimension vcrticale vpix. Pendant 
Tacquisition, les images sont transferees automatiquement par la carte d'acquisition dans le tableau qui 
leur est reserve en memoirc centrale de I'ordinateur. 
5 (801): le programme calcule les deux representations frequentiellcs SO et SJ qui different par 

Petat du commutateur. SO correspond au commutateur ouvert et SJ au commutateur fermc. 

Sty J] = ±{2I[0,0JJ] - /[1J0 J J] - /[2,0,/jD + - I[2fijj]) 

si[/,/] = ±(240.1,/./] - i[ujj] - i[2AJj]) + 7~(i[iXij] - /[2,u,/]) 

(8()2):Le programme calcule la valeur maximale mod^max = ^^^l^'Al atteintc P ar 

10 le module des elements du tableau SO. 

(803): le programme calcule 1c rapport moyen entre les deux representations frequentiellcs. II 
initialise rapport et nombrejvaleurs a 0 puis parcourt Tcnsemble des indices (ij) en testant la condition 
(804): |S0[/, j J £ 0,5 mod_ max 
Lorsque la condition est verifiee, il effectue (805): 

15 rapport* = ^--j 
nombrevaleurs + « 1 

Lorsque Tensemble des indices /\y a etc parcouru, le programme divise rapport par nombre jvaleurs 
(806) ce qui donne le rapport recherche. 

(807) le programme calcul la moyenne rapport _moy des valcurs de rapport obtcnues depuis 



/ 20 son lanccmcnt 

v 



(808) : le programme affiche les parties reelles et imaginaires ainsi que le module dc rapport et 
rapport moy. 

(809) Le programme reitere continument cette procedure pour permettre le reglage en continu. 
Le programme se tcrminc sur instruction de Toperateur. 

25 La position angulaire du rotateur de polarisation doit etre reglee pour que le module de rapport 

soil a pcu pres egal a 8. Lc programme est alors arrete et relance\ et au bout d'un nombre d'itcration 
suiTisant la valeur moyenne complexe rapport moy est nolec et servira de base dans la suite des 
operations. 

30 obtention de la constante K et reglage du diaphrag me ( 1 14) et du miroir (1 16) 

A' est la valeur maximale en pixels correspondant a la frequence spatiale maximale de Fondc 
j ou A est la longueur d'onde dans le milieu observe. II y a K pixels entre les frequences 0 el - Lc 
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pas en frequence suivant un axe est done -J- . Les frequences varient au total de - - a — par pas dc 



1 i r l _1 ^— *4 (Atnt Aa l_ O 

KX 

1 

KX 



Si N est le nombre total de pixels suivant chaque axe pris en compte pour la transformee dc 

.. a N * N 

5 Fourier, W valeurs de frequences sont prises en comptes. allant de - -^j a . 

Apres transformation on obtient N valeurs de position avee un pas en position egal a la moilie 
de l'inverse de la frequence maximale avant transformation. 

1 KX 

Le pas en position est done — = — 

1 2KX 

Si on considere deux points entre lesquels la distance en pixels est et la distance reclic est 

10 D„ e , . on a done: D wl = — d'ou K = -JIT^ La longueur d'onde a consider ici est la 

longueur d'onde dans 1c materiau soit si n est 1'indice du materiau et X v la longueur d'onde du laser 

dans le vide: X = — . On a finalemcnt: 
n 

K-JLJL D 

K -KD po D "*- 

Pour obtenir la constante K, on realise Timage d'un micrometre objectif pour lequel les 
1 5 distances reelles sont connues, puis on applique la formulc ci-dessus. 

Ceci est realise en utilisant un programme de focalisation qui sera reutilise ulterieurcment a 
chaque fois qu'une focalisation sur un tkhantillon est n^cessaire avant la prise d' image 
tridimensionnelle. L/algorithme de cc programme se trouve sur la figure 14. Ses Stapes principales sont: 

(1000) : Lc programme acquiert une image par la procedure acquiert images de la Fig. 13. On 
20 obtient done un tableau d'entiers (type unsigned char pour 8 bits) /[p.cjj] ou Tindice p variant dc 0 a 

2 correspond a Tetat de phase, l'indice c variant de 0 a 1 correspond a Tetat du commutateur (0=ouvcrt 
l=ferme) et les indices r et j variant de 0 a hpix-i et de 0 a vpix-\ correspondent aux coordonnees du 
pixel. 

(1001) :Un tableau de booleens H[ij] est genere: il est initialise a 0, puis pour chaque pixel, la 
25 valeur maximale attcinte sur les trois images correspondant a la position ouverte du commutateur est 

calculee. Si cette vaicur est egalc a 255 (valeur la plus clcvee du numeriscur), les images correspondant 
a la position fermee du commutateur doivent etre utilisces dans le calcul dc la frequence associec au 
pixel et a ses 8 voisins immediats, et lc tableau est mis a 1 pour ccs 9 pixels. 

(1002) : la representation frequentielle S[iJ] de nombres complexes est generee: 
30 pour chaque pixel, la valeur est generee selon les equations suivantes: 

S[/j]-[±(2/[a^ 
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^(l[lM[ij]jj] " l[2M[ijV >/]) j(l + (rapport jnoy - l)lf[ij]) 

Si vaut 1, la valeur complexe ainsi obtcnue est done multipliee par le nombre complcxc 

rapport jnoy obtenu lors de 1'operation de calibration du commutateur pour donner la valeur finale dc 
Telement de tableau S[ij], afin de tenir compte du decalage de phase et de l'absorption induits par le 
commutateur. 

(1003) : Le programme iimite le tableau S a des dimensions de 512x512. Le programme 
cfiectue ensuite optionnellemcnt Tune ou F autre, ou aucune, des deux operations suivantes: 

- un moyennage sur une largeur de 2, qui ramene le tableau Saun tableau S' de dimensions 256x256 
avec S?[ij] = S[2i2 j] + S[2i-+ 1,2 j] + S[2i.2/+ I] + S[2/ + U2j + 1] . Ce moyennage, couple a une 
reduction de rouverture du diaphragme (114), permet de diminuer le diarnetrc dc la zone obscrvee et de 
reduire le temps de calcul. II est equivalent a un filtrage passe-bas suivi d'un sous-echantillonnagc. 

- une limitation de 1 ^ensemble des frequences observees a un carre de 256x256 pixels avec 

ST[ij] = S[l 28 + 1\128 + j] . Ccci permet de diminuer le temps de calcul au prix d'une reduction de la 
resolution. 

Toutefois dans le cas present le programme n'effectue aucune de ces deux operations. 

(1004) : Le programme effectue alors la transformee de Fourier inverse du tableau ainsi obtenu. 

(1005) : 11 affiche le resultat sur Tecran, en extrayant le module, la partie reelle ou la partie 
imaginaire.Dans le cas present il affichera le module. Quelle que soit la variable affichee, le tableau de 
nombre reels correspondant est d'abord normalise, soit par rapport a la valeur moyennc, soit par rapport 
a la valeur maximale. Le programme ecrit egalement sur disque le fichier de reels correspondant. 
Lorsque la partie reelle ou la partie imaginaire sont representees, il est essentiel que le point d" impact 
du faisceau direct soit bicn en (256,256) sur V image de taille 512 x 512, faute de quoi une modulation 
devient visible. Lorsque lc module est represente, le point d'impacl exact du faisceau direct n'influencc 
pas notablemcnt le resultat. 

(1006) :Le programme recommence alors Tacquisition d'une nouvelle image, operant done en 
continu. II s'arrete sur instruction de Toperateur. 

Pour obtenir r image du micrometre, Tobjectif est d'abord mis en position a peu pres focalisec. 
lc micrometre ayant ete introduit commc objet. Dans un premier temps Tondc de reference est 
supprimec, les diaphragmcs sont ouvcrts au maximum, el lc programme dc visualisation dirccte en 
temps reel de V image rccuc sur le capteur est lance, les filtres en (103) et le rolateur de polarisation 
etant regies pour laisser passer suffisamment d'intensite pour saturer assez largcment le capteur au 
point d'impact direct du faisceau. L'objet est alors deplace a dans le plan horizontal a l aide du 
positionneur correspondant jusqu 1 a ce que Fimage constitute de nombrcux points intenses alignes, 
caracteristique du micrometre, apparaisse. Le micrometre est alors correctemcnt positionne sous 
l'dbjeclif. 
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Les diaphragmes en (107) et (1 14) sont alors regies pour une ouverture cT environ 8 mm. Les 
fiUres en (103) sonl alors regies pour que le maximum d'intensite sur 1c CCD soil a un niveau 
d'environ le quart de la valeur maximale du digitaliseur, soit 256/4=64. L'onde de reference est 
reintroduite. Le programme de focalisation est alors lance. Le diaphragme (1 14) doit etre ajuste de 
maniere a etre clairement visible sur r image affichee, tout en etant aussi ouvert que possible. Si F image 
n'est pas bien centrec, on peut ameliorer le centrage soit en modifiant roricntation du miroir (116), 
auquel cas il peut etre necessaire de reajusicr I'orientation du miroir (121), soit en modifiant la position 
du diaphragme (1 14). Le diaphragme en (107) doit etre ajuste pour que la zone obscrvee apparaisse 
uniformcment eclaire. Le programme de focalisation est alors arrete, Tonde de reference est supprimee 
et Tintensite du faisccau est reajustee comme precedemment. L'onde de refdrence est alors reintroduite 
et le programme de focalisation relance. 

L'objectif de microscope est alors deplace par le dispositif de focalisation de maniere a obtenir 
une bonne image du micrometre. Pour faciliter la focalisation, on a interet a visualiser une partic du 
micrometre ou des trails de differentes longueurs sont presents. Ceci limite les Tausses idealisations' 
dues a des phenomenes d'interferences en avant du micrometre. Entre deux deplacements il est 
necessaire de lacher le dispositif de focalisation manuel pour obtenir une image non perturbee par les 
vibrations. Lorsque une bonne image a ete obtenue, le programme est arrete et V image obtcnue est 
utilisee pour obtenir la distance en nombre de pixels entre deux traits, la distance metrique entre ces 
traits etant connue. Si la distance entre deux graduations separees par D rce i micrometres est sur 
Pimage ainsi obtenue de D pix pixels, si 1'indice moyen dans Techantillon est n (en general, n =L5) ct 
si la longueur d'onde du laser dans le vide est A v ( X v =0,633 micrometres) et si le nombre de points de 
la transformee de Fourier est N (N=512) alors on a: K = D reel , ou bien entendu D reel el A v 

sont dans la meme unite. 

reglage des diaphragmes 

L'image iridimensionnellc que Ton calculera a un cote de taille 256 pixels, ce qui permel dc 
limitcr les tailles de fichiers et les temps de calcul. Le reglage du diaphragme consiste a reutiliscr 1c 
programme de focalisation, cette fois avec ['option de moyennage intermediaire, et regler le diaphragme 
(1 14) de maniere a cc que son image soit clairement visible tout en 6tant aussi grande que possible. Lc 
diaphragme (107) est alors regie pour qu'il soit legerement plus ouvert que le minimum permettant un 
cclairagc r^gulier de la partic observee de Techantillon. 

aiuslement de la position du condenseur, obtcntion des caracteristi ques du cercle 
Lors de la procedure d 1 acquisition, rextremite du vecteur frequence d'eclairage doit suivre un 
ecrele (504), ce qui equivaut a ce que le point d'impact du faisceau d'eciairage sur lc capteur suive un 
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cercle correspondant. On appelle ici point d 1 impact du faisceau d'eciairage ie point illumine sur lc 
capteur (118) lorsquc seul le faisceau d'eciairage est present (lame transparente), et qui reste le point 
d'intensite maximale lorsque un dchantilion susceptible d'etre utilise par ce microscope est present. 
Pour une position donnee du condenseur (1 1 1) et de I'objectif de microscope (113) , le rayon du ecrele 

5 ne peut pas en general etrc choisi librement: ceci est une consequence des fortes aberrations induites par 
le condenseur pour les rayons non paraxiaux. Par ailleurs, si la position du condenseur est mal reglee 
transversatement, le faisceau d'eciairage son dc la zone observee de l'objet et est intercepte par lc 
diaphragmc (1 14). et le cercle ne peut pas etre parcouru. Si clle est mal reglee dans le sens dc Faxe 
optiquc, les caracteristiques du cercle ne sont pas optimales. 11 est done necessaire d'ajuster la position 

10 du condenseur et de calculer les caracteristiques du cercle. C'est ce que permet lc programme dont 

Talgorithme est detaille sur la Fig. 15. Le principe de ce programme est de deplaccr le point eclaire par 
le faisceau direct (point d'impact du faisceau d'eciairage) le long de rayons concentriques partanl du 
centre du capteur et d'afficher les vaieurs d'eclairement en chaque point eclaire successivement par le 
faisceau direct ce qui forme a Tecran une figure en 'etoile' representee schematiquement sur la Fig. 16. 

1 5 La position du condenseur peut alors etrc reglee pour que cette etoile soit la plus reguliere et large 
possible, et le programme determine les caracteristiques des cercles (centre commun G et rayons 
respectifs R p et R q ) delimitant la partie exterieure de I'etoile, qui permettent d'obtenir le cercle sur 
lequel se deplacera le point d'impact direct du faisceau d'eciairage lors de la procedure d'acquisition. 
Avant de lancer ce programme, Fonde de reference doit etre supprimee. Le miroir (109) doit etre 

20 deplace de fapon a ce que le point d'impact direct du faisceau soit au centre de Fimage. Les filtres en 
(103) doivent etre ajustes pour eviter la saturation du capteur et devront etre reajustcs pendant le 
fonctionncment du programme pour eviter la saturation sur 1 'ensemble de F etoile. Les etapes 
principales de ce programme sont: 

(1 100) : le programme calcule les vaieurs du rayon, indice par k, et de Tangle, indice par /. puis 
25 les coordonnees du point objectif. Les coordonnecs x,y de ce point sont exprimees en pas dc moteur (el 

non en pixels), x rcpresente le nombre de pas efifectues par le moteur (213) depuis sa position initiale 
(point eclaire au centre de Timage). v represente le nombre de pas effectu^s par le moteur (208) depuis 
sa position initiale. Lc pas du deplacement le long d'un rayon, pas_rayon s doit etre regie pour que les 
points de Tetoile soicnt suffisamment proches les uns des autres, et rincrcment angulairc^v angle 

30 vaut ln . I max vaut par exemple 16, k max est ajuste pour que le point eclaire pour I' indice 
/ max 

kjriax sorte du domainc ou il est visible. 

(1 101) : Lc programme deplace lc moteur vers lc point objectif 

(1102) : le point qui avait ete affichc sur i'ecran a Tetape precedente pour les memes indices 
kj est efface. Les coordonnees de ce point etaient rangees dans le tableau Ken \/\k,l\ 

35 (i 103): une image est acquise depuis le capteur. II s'agit la d^une acquisition simple, la vaicur 

de chaque pixel etant lue directemcnt sur le capteur. 
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(1 104) : le maximum declairement sur F image acquise est calcule et enrcgistre, ses 
coordonnees en pixels sont cnregistrecs dans le tableau ^contenant des couples d'entiers, en V[k.l), et 
sa valeur est enregistree dans le tableau W contenant des reels, en W\k.t]. Le point correspondent au 
maximum d'eclaircment est affiche sur Tecran (Fetoile est formee progressivement par ces points). 

( 1 105) : Lc programme determine la valeur maximale 

vmax = max (W[P* ']) atteinte sur la partie k > k_min du rayon. La valeur kjnin est une 

k jnin<>p<k jnax 

constante choisie de manicre a ce que la partie centrale intense de 1'etoile ne soit pas prise en compte. 
Le programme determine ensuite les entiers: 
kp = max[p e[k_ min, k_ max\ W[p, /] £ vmax. vcoef] 

kq = minj p e [*_ min. k max\ W[/>, /] 2> vmax. vcoef} 
ou vcoef=0,5 par exemple. II enregistre alors les coordonnees des points corrcspondants dans des 
tableaux Pet Q. 
P[l] = V[k P J] 

%l] = V[kq,l] 

1 5 Le programme affiche les valeurs de / et de vmax dans une fenetre de texte. 

(1 106) : Le programme calcule une approximation du centre du cercle et la range dans C[i]. 
C\I\ est calcule comme V intersection de la mediatrice de P{f) et P[!+l_max/2] avec la mediatrice de 
P\M maxIA] et P\!\l_max!M_maxI2]. La Fig. 17 illustre le principe de calcul sur Texemplc /_wox=8, 
pour le point C[0]. Le point C[0] est V intersection de deux droites dont Tune est la mediatrice du 

20 segment P[0] P[4] et r autre est la mediatrice du segment P[2)P[6). 

(11.07): Le programme calcule le centre G du cercle comme barycentre des points C[/J. Les 
coordonnees de G seront notees G x , G y 

(1108) : II calcule le rayon R p du cercle extcrieur comme valeur moyennc de la distance cntrc 
G ct les points P\l\ . et lc rayon R q du cercle interieur comme valeur moyenne dc la distance cntre G et 

25 les points Q\l\. 

(1 109) : II affiche G x , G y . R p , R q dans la fenetre de textc. 

(11 10) : II recommence Tensemblc de 1'operation. Le programme s^rrete sur demande de 
lopcrateur. 

Pour cette operation de reglage. le micrometre doit etre legerement deplacc de manicre a etre 
hors du champ d^observation, Tobjet etant alors equivalent a une lame transparentc. L'objectif doit etre 
maintenu en position focalisee. Dans un premier temps, le programme est utilise pour affiner le 
positionncment du condenseur. Celui-ci doit etre positionne transversalcment de manicre a ce que les 
points affiches forment une 'etoile' dont les branches ont toutes une largeur (notee L sur la Fig. 16) et 
une iuminosite proche, c ? est-a-dire que les valeurs des maximums affiches dans Tetape ( 1 105) doivent 
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etre proches pour toutes les branches. 11 doit etre positionnc dans la direction de laxe optiquc de 
manicre a ce que ie rayon R p du cercle cxterne soit aussi elevc que possible ct a ce que sur chaque 
branchc dc 1 'etoile, la zone cxterieure correspondant a des intensitds devees garde une largeur 
suffisante (par exemple supcrieure a 10 pixels), notee L sur la Fig. 16. Plus cettc /.one est etroite plus le 

5 programme d'acquishion d'image, qui utilise les coordonnees obtcnues pour lc cercle, risque de 

diverger en essayant de parcourir ce cercle. Pour cette premiere etape, il suffit d'assez peu d'angles <p . 
par exemple 8. Eventucllement P ouverture du diaphragme (107) peut etre augmentee, ce qui permet 
d'augmenter la largeur Dans un second temps, le nombre d'angles pris en compte est plus elevc afin 
d'obtenir des valcurs precises de G x ,G r R p , R q . Lc point G est le centre du cercle qui sera 

10 parcouru lors de Tetape d'acquisiuon par le point d'impact direct du faisceau d'eclairage. Le rayon R dc 
ce cercle est determine a partir des rayons R p et R q . Si A= R p - R q est faiblc (inferieure a 20 pixels 

R + R 

par exemple), on prend R = p .^ q . Si L est plus elevde, on peut prendre une valeur legerement 
inferieure a R p , par exemple R = R p - 10 . Si le condenseur est utilisd a une ouverture suiTisamment 
inferieure k son ouverture maximale, la partie externe de l'etoile rejoint le centre, la largeur L va done 
15 jusqu'au centre, et on peut la aussi utiiiser R = R p - 10 . 

enregistrement de Tonde de reference 

La connaissance de Tonde de reference est indispensable au calcul precis des valeurs 
complexes de Tonde parvenant au capteur. Celled doit done etre enregistnSe independamment dc la 

20 valeur constante de bruit moyen qui caracterise chaque pixel. A cet effet un programme specifiquc est 
utilise. Dans un premier temps, les ondes d'eclairage et dc reference sont supprimees et le programme 
enregistre T image de 'noir optique' qui en resulte sur le capteur CCD. II moyenne l'intensite obtenuc 
sur 100 images pour avoir un noir optique debruite. Dans un deuxieme temps Tonde de reference est 
retablie ct Tondc d'cclairage reste supprimee. Le programme enregistre V image resultante, la moyenne 

25 sur 100 images pour debruiter. Puis lc programme calcule la difference entre Timage dc Tonde de 

reference seuie et Timage de noir optiquc, ct enregistre Timagc rcsultante dans un tableau Jref\iJ] ou / 
varie de 0 a hpix-i et j varie de 0 a vpix-l. 

Focalisation sur Tobiet etudid 
30 Cette etape doit etre reiterce pour chaque echantillon dont on souhaitc obtenir une image. 

Londe d'cclairage est retablie. L'echantillon a etudier est mis en place. Le miroir (109) est regie pour 
que le point d'impact direct du faisceau d'eclairage soit au centre du capteur. Les filtres en (103) sont 
regtes pour qu'en Tabscence d'onde de reference Pintensite maximale recpue sur le capteur CCD soit 
d'environ 64. L^onde de reference est alors retablie. Le programme de focalisation est lance avec 
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roption de moyennagc intermediaire, et la position de I'objectif est ajustee a Taide du dispositif de 
focalisation pour obtenir une image nelte de ia zone d'interet dc rechantillon. 

Aiustement du ravon R 

5 Apres la phase de focalisation la phase d'ajustemcnt du condenseur doit ctrc repetee, sans 

modifier la position de P objectif et en laissant en place l'echantillon, et en utilisant la position fermec 
du commutateur optiquc. La position du condenseur n'est pas necessairement modifiec mais ia nouvclle 
valeur de R est notee et les filtres en ( 103) doivent etre ajustee de maniere a ce que rintensitc maximale 
rc<?ue sur chaque rayon (afTichee dans Tetape (1 104) de Talgorithme du programme de focalisation) soil 

10 d'environ le quart de la valeur maximale du digitaliseur, soit 64. L'onde de reference, qui doit ctrc 
supprimee pendant cetle etape, est rdtablie pour r etape d' acquisition. 

etape d* acquisition 

Gette etape permet d'acquerir les representations frequentielles bidimensionnclles a partir 
15 desquelles sera calculee la representation tridimensionnelle. Le point d 1 impact direct du faisceau 

d'eclairage est deplace suivant un cercle de centre G et de rayon R determines dans T&ape de reglage 

qui precede. Le pas angulaire pas est de ou nhim est le nombre de representations frequentielles 

bidimensionnelles a acquerir et a pour valeur la partie entiere de 2nR A chaque etape, lc programme 
calcule et enregistre dans un fichiQifich acquis une representation fnkjuentielle qui sera utilisce dans la 
20 phase de calcul tridimensionnel. La taille de l'image acquise depuis la camera est hpix x vpix, mais 
cette taille est divisee par deux pour obtenir des representations frequentielles de dimension he/ x veL 
selon lc principe de moyennage intermedial re deja utilise dans le programme de focalisation. 

Lc programme dacquisition est detaillee par Palgorithme de la Fig. 18 dont les etapes sont les 
suivamcs: 

25 ( 1 300): Les donnees de base de r acquisition sont enregistrdes dans le fichitr fich acquis: 

nbim: nombre de representations frequentielles bidimensionnelles 
hel=hpix 12 : nombre final de points dans le sens horizontal 
ve/^vpix/2: nombre final de points dans le sens vertical 

(1301) : temps d'altcnte permettant Tabsorption des vibrations creees par le mouvement du 
30 positionncur (110). Un temps d attcnte de 2s environ pcut convenir. 

(1302) : Le programme acquiert les images par la procedure acquiertjmages dc la Fig. 13. On 
obtient done un tableau dentiers (type unsigned char pour 8 bits) l\p.cjj] ou Tindicc p variant de 0 
a 2 correspond a retat de phase, Tindice c variant de 0 a 1 correspond a retat du commutateur 
(0=ouvert, l=ferme) et les indices i et j variant respectivement de 0 a hpix-} ct de 0 a vpix-l 

35 correspondent aux coordonnees du pixel. 
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(13()3):Un tableau de booleens H[ij] est gendre: il est initialise a 0, puis pour chaque pixel, la 
valeur maximale atteinte sur lcs trois images correspondant a la position ouverte du commutateur est 
calcul&s. Si cette valeur est egalc a 255 (valeur la plus dlevee du numeriseur), les images correspondant 
a la position fermee du commutateur doivent etrc utilisees dans le calcul de la frequence associee au 
pixel et a ses 8 voisins immcdiats. et le tableau // est mis a 1 pour ces 9 pixels. 

(1304): la representation frequentielle S[iJ] de nombres complexes est generec: 
pour chaque pixel la valeur est gen^rce selon les equations suivantes: 



i 



=(2/[o, H\ij]jj)- i[\,h[u]jj]- /[2,//[i\y],/,y]) 



+7 ) rU[UH[ij}jj]-l[2M[ij]jj]l [l+(rapport_mqy-\)H[u]) 

2^3Iref[i 9 j] y J 

10 Si vaut 1, la valeur complexe ainsi obtenue est done multipliee par le nombrc compiexe 

rapportjnoy obtenu lors de F operation de calibration du commutateur pour donner la valeur finale dc 
r element de tableau S[ij]. afin de lenir compte du decalage dc phase et de rabsorption induites par le 
commutateur. 

(1305) -Le point de module maximal du tableau S est determine. Ses coordonnees (i_max, 
15 jjnax) sont enregistrees. 

(1306) La position 'objectif est caiculee, Tangle valant k.pas: 
i_obj = G x + Rcod(k.pas) 

j_obj = G y + Rsiv(k.pas) 

ou ijobj et j__obj sont les coordonnees de la position objectif en pixels. 
20 Le deplaccmcnl des moteurs est calculc: 

pas /=(/ obj-i jnax).pas par _pixe! 
pas y obj-J max). pas par pixel 
ou pas j>ar j>ixel est le nombre de pas des moteurs par pixel de deplacement, determine 
cxpcrimentalement pendant lcs reglages. 
25 (1 307)- Un des moteurs permet de deplacer le point d'impacl direct du faisceau dans la 

direction dc Taxe /. II doit effectuer un nombre de pas pas J dans la direction correspondant a un indice 
i croissant pour la position du point d'impact direct du faisceau d'eclairage (pour pas J<0. il doit 
cfTcctuer un nombre de pas -pas i dans la direction opposee). L'autre moteur permet de deplacer le 
point d'impact direct du faisceau dans la direction de I'axc j. De meme, il doit effectuer un nombrc dc 
30 pas pas J. 

(1308)- Un tableau de frequences moyennees est genere: 
Chaque dimension du tableau initial S est divisee par deux pour donner un tableau M k avec 
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(13()9)-Le point de module maximal du tableau M k est determine, ses coordonnees imax k . 
ymar* et sa valeur/wajr_ ™oy = Kf k \tmax k ,jmax k ] = ^^^(^ * [M) sont enregistrees. 

(1310) -Les elements du tableau M k sont normalises cn les divisanl par maxjnoy. 

u \ 1 M M 
1 J wax_ m«v 

(1311) - imax k % jmax k et la representation frequentielle M k sont enregistres dans lc fichier 
fichjacquis. 

(13 12) - L'algorithme se termine quand Tangle vaut In . Le moteur revient alors a sa position 
initiale et le fichier d'acquisition fich acquis est ferm£. 



10 



calcul tridimensionnel 

La procedure d'acquisition a gendre des representations frequentielles bidimensionnelles M k s 
ou Tindice k represente le numero d'ordre de chaque representation . Ces representations ayant ete 
moyennees sur une largeur de 2. les valeurs de C x X y et K doivent etre divisees par 2 pour 
15 correspondre au nouveau systeme de coordonnees. Ce sont ces valeurs divisees qui sont utilisees dans la 
suite. L/ensemble des representations bidimensionnelles peut etre considerc une representation a trois 
dimensions (i.j,k) dans lequel Tindice k represente Tindice image et les indices / ety represented les 
coordonnees cartdsiennes de chaque representation. Dans la k -ieme representation frequentielle 
bidimensionnelle: 

20 la relation entrc les coordonnees ct la frequence caracteristique est: 

La relation entrc les coordonnees et la frequence d'eclairage est: 

KJfe = {ima*k -C x Jmax k -C y ^K 2 -(imax k -C x ) 2 -{jmax k -C y ) 2 j 
La relation entre les coordonnees cl la frequence totalc est done: 



25 KVi = - imax k , j - jmax k . ^K 2 -(/ " C x f - (j -C y ) 2 - ^K 2 -{imax k -C x ) 2 -{jmax k -C y ) 2 j 

La procedure de calcul tridimensionnel consiste en principe a generer dans un premier temps une 
representation tridimensionnelle sous forme dun tableau Fde dimensions fdim xfdim xfdim avec 
fdim=256 t puis a cn effectuer la transfonnee de Fourier pour obtenir un tableau U de mcmes 
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dimensions correspondent a la representation w(r), ct dans lequel les indices correspondent done au 
rayon-vecteur. La representation F est un tableau dont les indices representent les coordonnees de 
KAf t , le zero etani ramene par translation au point de coordonnees (fdimlljtiimlljdimil). On obtient 
done, a partir de cliaque point dc coordonnees ij d'une representation bidimensionnellc M k , un point 
de la representation tridimensionnelle par: 
F[ni.>y\nk] = M k [ij] 
avec: 

(dim 

ni - i - imax k + — — 
fdim 

nj = j-jmax k + —— 



nk = ^K 2 ~{i-C x ) 2 -{j~C y ) 2 -^K 1 -(imax k -C x ) 2 -(jmax k -C y ) 2 +^ 



2~ 



Lorsque un point du tableau F, de coordonnees (ni,nj,nk) y est obtenu succcssivement a partir dc 
plusieurs representations bidimensionnclles distinctes, la valeur de F retenue en ce point est la moyenne 
des valeurs obtenues a partir de chacune de ces representations bidimensionnelles. Lorsque ce point 
n'est jamais obtenu, on y affecte une valeur nulle. On peut egalement retenir comme valeur de F la 

15 somme des valeurs obtenues: ceci equivaut a I'approche citee dans la description de l'invention et 

consistant a moyenner sur I 'ensemble des representations bidimensionnelles y compris celles ayant une 
valeur nulle au point considere. Ce n'est toutefois pas cette solution qui est retenue ici. 

Lorsque le tableau Fa ete genere, le tableau U peut etre obtenu par transformee de Fourier 
tridimensionnelle inverse. Pour des raisons pratiques, une partie dc la transformation dc Fourier sera 

20 effectuee au fur et a mesure de la generation du tableau F. qui ne sera done jamais veritablement 
genere. 

L'algorithme du programme de calcul est decrit Fig. 19. Le programme opere en cinq ctapes. 
Chaquc etape utilise un fichicr d'entree stocke sur le disque dur de Tordinateur et genere un fichier de 
sortie egalement stocke sur lc disque dur, dont le nom est en italique sur la figure. Le programme 

25 pourrait thcoriquement effectuer de maniere plus directe les operations nccessaires a la generation dc 
Pimage tridimensionnelle, mais la taille des fichiers impliques est trop elevee pour qu'ils puisscnt etre 
entierement contcnus dans la memoire d'un ordinateur. II est done necessaire de gerer leur stockagc sur 
disque. La lecturc/ecriture d 'elements de fichiers sur disque sc faisant plus rapidement si les elements 
sont stockes de maniere contigue, le programme doit etre concu pour nc lire ct ecrire sur lc disque que 

30 des blocs d'unc taille suffisante. C'est cc que permet 1 'algorithme decriL alors qu'une methodc directe 
necessiterait des lectures/ecritures a des addresses non contigues et ne serait pas praticable en raison du 
temps perdu en acces disque. Les etapes de Falgorithmc sont les suivantes: 

(1400)-centrage des representations frequentielles bidimensionnelles: 
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Cette procedure consiste a translater les representations bidimensionneiles pour passer (Tunc 
representation dans le systeme de coordonnees (fj\k) a une representation dans le systeme (ni.nj.k) avcc 
ni-MnuBtfdimn.rtf-Jrfmax+fdim/1. les dimensions suivant ni ct nj du tableau ainsi gcner£ etant dc 
fdim xfdim. L'algorithme de cette panic est represent figure 21. Les etapes sont les suivantes: 

(1600) : les valeurs de he!,vel,nbim sont lues dans le fichier ft ch_acquis. 

(1601) : la representation ftequentielle M k correspondant a Tindice k est transferee en memoirc 
centrale avec les valeurs correspondantes imax k ti jmax k 

(1602) : Une representation frequentielle translate T k de dimensions/^//* xfdim est generee 



avec: 



10 



. T k [ni, nj] = M k [ni - fdim 1 2 + imax k , nj - /dim / 2 + jmax k ] 

lorsque 0 <£ /i/ - fdim 1 2 + imax k £ hel - 1 et 0 £ nj - fdim 1 2 + jmax k <Lhel-\ 



( . T k [ni.nj] = 0 dans les autres cas 

La Fig. (20) montre une image d'origine avec son point de coordonndes ( imox k , jmox k ) et un 
dessin arbitraire autour de ce point, et la nouvelle image translatee. 
15 (1603): Les valeurs de imax k , jmax k , et la representation T k sont enrcgistrees dans le fichier 

flehjeentre dans Tordre suivant: 

imax k Jmax k ,T k [0,0], T k [1,0], T k [fdim-l,0], 

imax k . jmax k . T k [047-414 T k [fdim- 1,1], 



20 imax k , jmax k , T k [0,fdim)T k [l, /dim], T k [fdim - 1, fdim- 1] , 

Notons que /max* et /mar* sont repetes a chaque ligne de T k dans le fichier fich_cenire. 
Ceci permct que ces donnees, indispensables au changement d'indice k-+nk, restent disponiblc apres 
Toperation d'echange des axes. 

(1604): le processus est recommence tant que k est inferieur a nbim. 
25 (140U-Echange des axes olan^ et Migne' 

Le changement d 1 indices / -* ni et j -> nj ayant ete efTectue, il restc a faire le changement d'indicc 
k : -+ nk . Pour effectuer ce changement d'indice cn un temps raisonnable il est necessairc de pouvoir 
charger erememoire centrale, rapidemcnt, un plan (ni,k). Pour que cette operation soit possible, il faut 
prealablcment echangcr les axes k ct nj. Ccsl ce que fait cct algorithme. 
30 Le fichiery/cA centre prccedemmcnt cree est relu. Tordre des donnees est modifie et un fichier 
fich ech! est ecrit dans lequel les donnees sont dans Tordre suivant: 

» imax 0 , jmax Q , T 0 [0,0], T 0 [ 1.0], 7b [fdim - 1,0], 

[0,0], 7 1 , [1,0], 7",[^*m-l,0], 
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imax nbim _i Jmax nbim _ x , r„ Wm _,[0.0]. T„ 6/m _, [1,0], T nbim ^ [fdim- 1,0], 

imax 0 Jmax 0 , 7' 0 [0,l], 7* 0 [U], 7b[yHwt - 1,1], 

imaxi , jmax x , 7\ [0,1], 7, [l,l], 7', [/<A/w - 1.1], 

5 , jmax nblm _ | . 7„ Al ,„_, [0.1], 7*„{, lm -i [1,1], 7 , „ 6lm _ l [fdim - 1,1], 



i>»ar 0 , jmax 0 , T 0 [0, - l], T 0 [1, /dim - l] T 0 [fdim - 1, fdim - l], 

//nor, , jmaxi . 7\ [0, /<0m - l], 7\ [ 1, fdim - 1] 7*, [/rflm - 1, fdim - l], 

imax nbim ^Jmax nhm . x , ^^[0, ^*m- 1] F 7^,(1, >»»-l] 7^m-l !./*»- *] 

Cest-a-dire que les axes wy et Jfc sont echanges. Cctte operation d*echange des axes est eflectuee bloc par 
bloc. La Fig. (22) represente symboliquement (1700) le contenu d'un fichier tridimensional 
correspondant a des indices i j'X range en memoire plan horizontal par plan horizontal et chaque plan 
etant range ligne par ligne, une ligne etant dans le sens de la profondeur sur le dessin. Les axes ij\k et 
Torigine O du repere sont precises sur le dessin. Le contenu du fichier obtenu par inversion des axesy et 
k est represente en (1702). Le transfert des donnees d'un fichier a l'autre se fait bloc par bloc, lc bloc 
(1701) etant copie en memoire centrale puis transfere en (1703). La lecture et recriture des blocs se fait 
plan horizontal par plan horizontal, un plan horizontal dans le fichier lu ne correspondant pas a un plan 
horizontal dans lc fichier ecrit. La taillc du bloc est la taille maximum qui puisse tcnir dans la memoire 
interne disponible dc Tordinateur. L'opcration est repetee pour tous les blocs, les blocs situes *au bord' 
n'ayant pas en general la meme dimension que les autres. Cctle maniere d'operer pcrmet d'eflectuer 
Fechange des axes avec un ordinateur dont la memoire interne a une taille inferieurc a celle des fichicrs 
utilises. 

La procedure en language C (Microsoft C/C++ sous Windows 95) qui permet cctte operation est la 
suivante: 

void echangc_axes(FILE* rcad_file ? FILE* write_file,int ktot.int jtot,im itotjnt memoryjimit ) 

{ 

30 int knumJnum,bknum ? bjnum.kefTjerT.k j,bk,bj; 
char* buff: 

knum=(int)sqrt(((double)memor>'_limit)/((double)itot)); 
jnum=knum; 
35 bufr^(char*)malloc(itot*knum*jnum); 
bknum=ktot/knum; 
if (knum*bknum!=ktot) bknum+=l; 
bjnurn=jtot/jnum; 
if (jnum*bjnum!=jtol) bjnum+=l; 
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for (bk=0;bk<=bknum- 1 ;bk++) 
for (bj=0;bj<=bjnum-l;bj++) 

5 * jf (bk=(bknum-l )) kerT=ktot-knum*(bknum- 1 ); else keff=knum; 

if (bj==(bjnum-l)) jcfT=jtot-jnum*(bjnum-l); else jeff=jnum; 

for (k=0;k<=keff- 1 ;k++) 
for 0=Ou < =jeff-lU ++ ) 
10 { 

fscek(readJTile,(bk*knum+W 
freadCbufiF+k+jefi^itot+j+itot, L itot, read_filc); 

} 

15 for(j=0J<=jeff-lJ+4-) 

for (k=0.k<=keff-l;k++) 

( f se ek(writejlle,(bj^num 

fwrite(buff+k*jefT*itot+j*itot 1, itot, write_file); 

20 } 
} 

free(buff): 
> 



25 Les parametres a passer lors de Tappel de la procedure sont: 
read_file: pointeur sur le fichier fich centre 
write file: pointeur sur le fichier fichjechl 
ktot: nombre total d' images nbim 

jtot: nombre total de lignes dans une representation frequentielle :_/<///?* 
30 itot: taille en octets d'une ligne: /^/m*sizeof(complexe)+2*sizeof(int), ou sizeof(complexe) dcsigne la 
taille en octets dun nombre compiexe ( T k [ij] par exemple) et sizeof(int) la taille en octets d'un 
nombre entier ( imax * par exemple). 

memory Jimit: taille maximalc en octets de la memoire a acces aleatoire (RAM) dont dispose la 
procedure pour stocker les blocs. 
35 Les fichiers doivent etre ouverts en mode 'commited^ e'est-a-dire que les operations de 

lecture/ccriture sont eflectuecs directement depuis ou vers le disque dur, sans bufferisation intcrmediairc 
en memoire centralc. 

(1402V phase dc calcul. 

Cctte phase de calcul a pour objet de remplacer Tindice ' image * k par Tindicc nk donne par laformule: 

40 nk-£j± + jK 2 -(i-C x ) 2 -{j-C y ) 2 - jK 2 -(lmax k -C x ) 2 -(jmax k -C y ) 2 

ou ij sont les coordonnees dans le repere dorigine (avant centrage). Dans le repere centre on a done: 
nk = ^K 2 - [ni - fdim / 2+ imax k ~C X ) 2 -{nj - fdim / 2 + jmax k -C y ) 2 
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-t)K 2 ~(imax k -C x ) -(jmcx k -L y ) + — 

Lorsque les memes indices (ij\k) sont obtenus plusieurs fois par remplacement de 1' indice image, la 
valeur retenue pour Felement correspondant de la representation frequentiellc tridimensionnellc est la 
moyenne des valeurs pour lesquelles les indices (ij,k) sont obtenus. 

5 Un point de la representation frequenticlle tridimensionnelle, de coordonnees (ni>nj,nk) , ne 

peut etre obtenu par ce changement de coordonnees qu'a partir d'un plan (jii.fc) donne et correspondant 
a son indice nj. Les plans (ni,k) peuvent done etre traites independamment les uns des autres. Lorsque 
dans un plan (wi.it) on a remplace I'indice k par t'indice nK il est done possible d'cfTectuer directcment 
la transformee de Fourier bidimensionnelle inverse de ce plan avant de passer au plan suivant. C'cst ce 

10 que fait cette partie du programme, dont ralgorithme est detaille sur la. Fig. 23. Ses dtapes principals 
sont: 

(1800) les elements suivants sont lus dans le fichior fichechl et transferee en memoire 

interne: 

imaxQ , jmax 0 , F 0 [0, nj\ T 0 [lnj] y T 0 [fdim - 1, nj\ 

15 imax\Jmax x , 

imax nhim - X Jmcx nhim - X . T nbtm ^ x [0, nj\ T nbim ^[lnj] T nbim _ x [fdim - Inj] 

(1801) les tableaux D nJ et Poids. de dimensions fdim xfdim, sont initialises a 0. 

(1802) La condition suivante est testee: 

20 (ni-fdim/2+imax k -C x ) 2 + (nj - fdim/ 2+ jmax k -C y ) 2 <[ K ~j 

ou o est Touverture du microscope, legerement surevaluee par precaution, ct n V indice de 

1 2^ 

rhuile optique et de la lamelle utilises, soit a peu pres. ^ = j— . Si la condition est vraie, le point 

correspond a un vecteur frequence ne sortant pas de Pouverture dc Fobjectif et est done dans la zone 
observable. 

25 (1803) La valeur nk est calculee par la formule: 

nk = ^K 2 - (ni - fdim / 2 + imax k -C x ) 2 -{nj - fdim / 2 + jmax k -C y ) 2 

-^K 2 ~{imax k -C x ) 2 ~{jmax k -C y f + fdim 1 2 

( 1 804) La valeur correspondantc de frequence est additionnee au tableau D nj L'element 
correspondant du tableaux des poids, qui sera utilised pour calculer la valeur moyenne, est increment^ 

30 D„jlni,nk)+ = T k [ni^nj] 
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Poicblni,nk]+ = 1 

(1805) Lorsque r ensemble des indices ni.nk a ainsi etc parcouru, le programme parcourt 
Pensemble des indices ni et nk en testant la condition Poids{m,nk) *0 et a chaque fois que cette 
condition est realisee il effectue: 

D nj [ni,nk] 

5 ' D„Ani*nk] = — = -r 

njl 1 Poid^ni.nk] 

(1806) le programme effectue une transformee de Fourier bidimensionnetle inverse du tableau 

D nj . 

(1807) il stocke le tableau transform* dans le fichier de sortie fich caic dans Tordre suivant: 
^njM,D nJ [lfi] D nJ [fdim-l,0l 

10 D n/ [04D^[l4:....D„ y [/^-l,l], 



25 



D nj [0, fdim - 1], D nj [1, fdim - l], D„, [/dim - 1, /tfni - i] 

(1403V2eme echange des axes: 
II reste a ce niveau a effectuer une transformee de Fourier inverse a une dimension suivant I'axe nj. 
1 5 Afin de pouvoir effectuer cette transformation en un temps raisonnable, il est necessaire d'echanger 
prealablement les axes nj et nk. Le programme pourra ensuite charger en memoire centrale des plans 
(ni, ny) complets pour les traiter. 

Dans le fichier fich calc les donnees sont rangees dans Tordre suivant: 

A) [0,0], D 0 [1.0], D 0 [fdim - 1.0], 

20 Z) 0 [OJ],D 0 [1,1] D Q [fdim-U\ 



D 0 [0, fdim - 1], D 0 [ 1, fdim - l] D 0 [fdim - 1, fdim - l], 

Dj [0,0], D x [1,0], D x [fdim - 1,0], 

D^O,!], D,[U], D x [fdim-\,\\ 



Dj [0, - 1], D\ [ 1, w - 1], D x [fdim - 1 fdim - i], 



£ [0,0], D ya /m _i [ 1,0] D fdim ^[fdim - 1,0], 

30 D^ /m _, [0,1], Dya^, [1.1]. D^., [/<*« - 1,1], 



Dfifo-x [0, JHfm - 1], D fdtm _ x [ 1, /rf#m - 1], D fdim ^ [fdim - 1, fdim - l] 
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Ce fichier est relu. et un fichier //c/? ech2 est genere, dans lequel les donndes sont reecritcs dans Fordre 



A,[0,0] ? £b[l,0], Do[fdim - 1,0], 

D,[0,0], D { [1,0], D x [fdim- L0], 

u fdxm-\ 

[0,0], D^JKO], D fdim - X [fdim - 1,0], . 

Bb[ai], A>[U], A>[.A*« - U]. 

D, [0,1], Z?i [1,1], D, [y»m - 1,1], 

0>Wm-l [04], D/Wm-j [U] D>Wm-l[>a« " l 4 



Z) o [0, /rfim - 1], D 0 [l, yam - 1], DolyWiw " *■ ^ " ! ]> 

£>, [0, yam - 1], D, [1, /dfm - 1], Di [fdim - 1, yam - l], 

Dfitim-l [0,/dim - 1], Dfd im -\ [l, ^tfm- l], DfiHm-\ [fdim - 1,/rf/m - 1] 

Cet echange des axes nj et /?£ est effectue par blocs comme le precedent echange d'axes. La memc 

procedure est utilisee, les parametres a passer &ant: 

read file: pointeur sur le fichier ftchcalc 

write_fiie: pointeur sur le fichier fichj>ch2 

ktot: fdim 

)tot: fdim 

itot: taillc en octets d'une ligne./rf//w*si/xof(cornplexe) ou sizcof(complexe) designe la laillc cn octets 
d'un nombre complexc. 

memorylimit: comme precedemment, la taille en octets de la memoire disponible. 

f 1404)-Derniere transformation de Fourier 
Cettc procedure consisted effectucr la transformation de Fourier inverse suivant Taxe nj. II s'agil d un 
traitement iteratif sur Tindice nk. L/algorithme de cette partie de programme sc trouve Fig.24. Ses 
etapes essentielles sont les suivantes: 

(1900): 1c programme charge en mdmoire interne les valeurs: 
D 0 [0.nk\ D Q [lnk\ D 0 [fdim - l.nk], 

D\ [0, nk], D x [lnkl D x [fdim - 1, nk\ 



suivant: 
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D fdm - x [0,nk\ Dfium-i[l»*] D fdim ^[fdim- Ink] 

(1901) :Le programme genere le tableau E nkjni a une dimension: 

II efifectue la transformee de Fourier inverse de ce tableau, generant le tableau E nk?nt 

(1902) : il enregisire les resultats dans le fichier fich_rep, dans Iordre suivant: 

£„*,Q [0], EnkA [0], Enkjdim- 1 [<>]> 



E nk%0 [fdim-\\ E nk +[fdim-\\....E r * Jdim ^ 
10 Le fichier ainsi genere contient alors la representation tridimensionnelle de l'objet au format lJ\nUnj,nk\ 
ou Pelement complexe U\ni t nj,nk\ est range dans le fichier fich rep a Taddressc 
{nk*fdim*fdim+nj*fdim+ni) comptee a partir du debut du fichier, Taddressage se faisant par elements 
de type 4 nombre complexe' . 

15 visualisation 

Le tableau U ayant ete genere, on peut en visualiser le contenu. La visualisation la plus simple 
consiste a en extraire une coupe, un des indices etant fixe a une valeur constants Sur cette coupe, on 
peut visualiser la panic reelle. la partie imaginaire ou le module. Un autre mode de visualisation 
consiste a en extraire des projections. Par exemple pour une projection de la partie reelle selon V axe nk. 

20 on representcra le tableau V avec 

V[m,nj] = ^[R^(u[ni f nj\nk])] 2 , ou Re(x) designe la partie reelle de x, la somme sur nk dtant prise 

entre deux plans ' limite' scion ce que Ton veut representer. 
Variantes : 

25 D'autres modes de realisation sont bien entendu possibles et la description ci-dessus nest pas 

limitative. Entre autrcs possibility on peut utiliscr des longueurs d'onde plus courtes, des captcurs plus 
rapides ou ayant un nombre de pixels different, on peut generer unc representation tridimensionnelle de 
taille diflcrente, on peut utiliscr un autre dispositifd' attenuation controlec, on peut utiliscr moins dc 
positionneurs et compter sur une realisation plus soignee, on peut modifier le systeme optique en 

30 adoptant unc configuration physique differente, par exemple en supprimant certains voire tous les 

elements optiqucs presents entre l'objet et le capteur. ou encore cn utilisant des ondes d'eclairagc non 
planes. 
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Possibtlites duplication industriellc: 

Ce microscope peut etre utilise au lieu dc microscopes a transmission usuels, dans le domainc 
de ia biologie ou de la metrologie. II a une meilleure resolution que Les systemes connus et permet unc 
visualisation dc Tindice local dans rechantillon. 
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REVENDICATIONS (1/2) 



10 



1 - Microscope calcuiant numcriquement une representation tridimensionnelle d'un objet observe a 
partir de plusieurs enregistrements de l'onde issue dudit objet et recue sur une surface de reception, 
different entre cux par l'onde eclairant ledit objet, caracterise par te fait qu'il comporte des moyens pour 
determiner une representation de chaque onde d'eciairage dans l'objet. 

2- Microscope selon la revendication 1, caracterise par le fait que ledit calcul numerique de ladite 
representation tridimensionnelle comporte une etape de calcul d'une representation normahsee de 
chaque onde issue de l'objet par rapport k l'onde d'eciairage correspondante. 

3- Microscope selon la revendication 2, caracterise par le fait que ledit calcul numerique de ladite 
representation tridimensionnelle comporte une etape de moyennage desdites representations 
normaiisees. 

^ 15 4-Microscope selon une des revendications 1 k 3 , caracterise par le fait que lesdites ondes d'eciairage 

sont des ondes planes different entre elles par leur direction. 

5-Microscope selon les revendications 2 et 4 caracterise par le fail que chacune desdites representations 
normaiisees est obtenue dans le domaine frequentiel en divisant V ensemble de la representation 
20 complexe de l'onde issue de l'objet par la valeur de l'onde d'eciairage sur sa-frequence unique, puis en 
translatant I'ensemble de ladite representation de maniere a ramener au zero la frequence unique de 
l'onde d'eciairage. 



25 



6-Microscope selon les revendications 4 et 5, caracterise par le fait que ladite representation 
tridimensionnelle de l'onde issue de V objet est obtenue par moyennage desdites representations 



normaiisees. 



7- Microscope selon la revendication 6, caracterise par le fait que ledit moyennage est effectue dans le 
domaine frequentiel et est limite en chaque point aux representations frequentielles normaiisees ayant 

30 une valeur non nulle au point considere. 

8- Microscope selon une des revendications 4 a 7, caracterise par le fait que Textremite du vecteur 
frequence de ladite onde d'eciairage parcourt un cercle. 

35 9-Microscope selon une des revendications 4 a 8, caracterise par le fait que les variations de direction dc 
ladite onde d'eciairage soni obtenues par reflection sur un miroir mobile renvoyant le faisceau 
d'eciairage vers un condenseur qui le redirige vers la zone observed de l'objet, un systeme optiquc 
permettant de focaliser ladite onde d'eciairage sur la surface focale objet dudit condenseur. 

40 10-Microscope selon unc des revendications 4 a 9, caracterise par le fait qu'il comporte un diaphragme 
permettant dc limiter la largeur du faisceau d'eciairage. 

1 1- Microscope scion une des revendications 1 a 10, caracterise par le fait qu'un dispositif d'attenuation 
controlec introduit sur le trajet du faisceau d'eciairage permel de faire varier Tintensite du faisceau 

45 d'eciairage et par le fait que chacune desdites representations normaiisees d'ondes issues de I'objet est 
obtenue k partir de plusieurs enregistrements de Tonde issue de Fobjet differant entre eux par l'intcnsite 
de l'onde d'eciairage. 

1 2- Microscope selon une des revendications 1 a 1 1, caracterise par le fait que les donnees 

50 tridimensionnclles sont stockees sur un support a acces sequentiel pendant le traitement et que le calcul 
de la representation tridimensionnelle est effectue plan horizontal par plan horizontal et comporte des 
phases d'echange d'axes ramenant provisoirement i'axe vertical dans un plan horizontal. 
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13-Representation tridimensionnelle. image bidimensionnelle ou image stereoscopiquc d'un objet 
observe, caracterisec en ce qu elle esi obtenue a l'aide d'un microscope selon une des revendicauons 1 a 

12. 



( 



( 
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valeur initiale de diff_haut,diff_bas 
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Acquisition des images 



602 
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calcul des representations frequentielle SO et SI 



recalagc cn phase de SI 



calcul de la valeur maximale du 
module sur SO 



ecart=Q* nombre valeurs=0J=0 
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non 



ecart+ = \SQ[ij]-S)[i,j]\ 
nombre vdleurs* = 1 
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ecart^ecart/nombre valeurs 
stockage d'ecart dans un tableau 
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611 




affichagc de la valeur de (dijfhaut.diffjbas) 
correspondant au minimum d' ecart 



> t 

Fig. 10 
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Regie tension du piezoelectrique a 0 



attend 20 ms 



Regie tension du piezoelectriqu e a Umax/2-d\ffJ)as 
3 



attend 40 ms 



met le processus en priorite maximale 
lance V acquisition 



T 



attend fin d" exposition lere image, 
stockee dans un tableau l\QS)J,j\ 

E 



Regie tension du piezoelectrique a Uma x/2-diffJ>as+d\ffJ)as/2 



attend fin d'exposition 2e image, 
stockee dans un tableau J\0,\JJ\ 



Regie tension du piezoelectrique a Umax/2 



attend fin d'exposition 3e image, 
stockee dans un tableau7[l,0,'\/J 



Regie tension du piezoelectrique a Umax/2+(diJjr_bas+diff_haut)/4 | 



attend fin d'exposition 4c image, 
stockee dans un tableau /[ 1 , 1 J J] 



X 



Regie tension du piezoelectrique a Umax/2 +diff_haut 



attend fin d'exposition 5e image, 
stockee dans un tableau J[2,Qj,j] 



Regie tension du piezoelectrique a Umax/2 \ dijf haut+diff hautll 




1 

attend fin d' exposition 6e image, 
stockee dans un tableau 7[2, t Jj\ 





regie la tension du piezoelectrique a Umax 



remct le processus en priorite normale 

1 ~ 



Fig. 1 1 
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Fig. 12 
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Regie tension du piezoelectrique a 0 



attend 20 ms 



Ouvre commutateur 
Regie tension du piezoelectrique a Umax/2-diJf_bas 



attend 40 ms 



met lc processus en priorite maximale 
lance V acquisition 



attend fin d'exposition lere image, 
stockcc dans un tableau /fOjO,/ J] 



ferme commutateur 



attend fin d' exposition 2c image, 
stockee dans un tableau I[0AJJ] 



Ouvre commutateur 
Regie tension du piezoelectrique a Umax/2 



attend fin deposition 3e image, 
stockee dans un tableau I\ 1 fijj] 



ferme commutateur 



attend fin d 1 exposition 4e image, 
stockee dans un tableau /[ 1 , 1 J J J 



Ouvre commutateur 
Regie tension du piezoelectrique a Umax/2 * cliff haut 



attend fin d'cxposition 5e image, 
stockcc dans un tableau /[2,0,'Vl 



fcrmc commutateur 



attend fin deposition 6c image, 
stockee dans un tableau I[2AJJ] 



regie la tension du piezoelectrique a Umax 



remet le processus en priorite normal e 



T 



Fig. 13 
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commutateur doit ctre utilisee 



1002 



1003 



1004 



calcul de la representation 
frequenuelle 



T 



limitation de taille et 
modifications optionnelles 

i 



Transformation de Fourier inverse 
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Fig. 14 
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acquiert image 
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Fig. 15 
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Fig. 18 
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fich centre 
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dc transformee de Fourier 
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fich rep 



Fig. 19 
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(1600) 
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calcul de la representation 
frequentielle centree 
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Fig.21 
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